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ABSTRACT 



'Water In Star-forming regions with Herschel' (WISH) is a key program on 
the Herschel Space Observatory designed to probe the physical and chemical 
structure of young stellar objects using water and related molecules and to follow 
the water abundance from collapsing clouds to planet-forming disks. About 80 
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sources are targeted covering a wide range of luminosities -from low (< 1 L0) 
to high (> 10^ Lq)- and a wide range of evolutionary stages -from cold pre- 
stellar cores to warm protostellar envelopes and outflows to disks around young 
stars. Both the HIFI and PACS instruments are used to observe a variety of 
lines of H2O, H2^^0 and chemically related species at the source position and in 
small maps around the protostars and selected outflow positions. In addition, 
high-frequency lines of CO, ^^CO and C^^O are obtained with Herschel and are 
complemented by ground-based observations of dust continuum, HDO, CO and 
its isotopologues and other molecules to ensure a self-consistent data set for 
analysis. An overview of the scientific motivation and observational strategy of 
the program is given together with the modeling approach and analysis tools 
that have been developed. Initial science results are presented. These include 
a lack of water in cold gas at abundances that are lower than most predictions, 
strong water emission from shocks in protostellar environments, the importance 
of UV radiation in heating the gas along outflow walls across the full range 
of luminosities, and surprisingly widespread detection of the chemically related 
hydrides OH"*" and H2O+ in outflows and foreground gas. Quantitative estimates 
of the energy budget indicate that II2O is generally not the dominant coolant 
in the warm dense gas associated with protostars. Very deep limits on the cold 
gaseous water reservoir in the outer regions of protoplanetary disks are obtained 
which have profound implications for our understanding of grain growth and 
mixing in disks. 

Subject headings: star formation, young stellar objects, protoplanetary disks, 
astrochemistry, molecular processes 



Introduction 



As interstellar clouds collapse to form new stars, part of the gas and dust are transported 
from the inf ailing env elope to the rotating disk out of which new planetary systems may 
form (jShu et al.l Il987| ). Wate r has a pivotal role in these protostellar an d protoplanetary 
environments (see reviews by ICernicharo &: Crovisierl l2005l iMelnickl 120091 ) . As a dominant 
form of oxygen, the most abundant element in the universe after H and He, it controls the 
chemistry of many other species, whether in gaseous or solid phase. It is a unique diagnostic 
of warm gas and energetic processes taking place during star formation. In cold regions, water 
is primarily in solid form, and its presence as an ice may help the coagulation process that 
ultimately produces planets. Asteroids and comets containing ice have likely delivered most 
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of the water to our oceans on Earth, where water is directly associated with the emergence 
of hfe. Water also contributes to the energy balance as a gas coolant, allowing clouds to 
collapse up to higher temperatures. The distribution of water vapor and ice during the entire 
star and planet formation process is therefore a fundamental question relevant to our own 
origins. 

The importance of water as a physical diagnostic stems from the or ders of magnitude 



variations in its gas phase abundance between warm and cold regions (e.g., ICernicharo et al. 



199fl, Ivan Dishoeck fc Helmichll 



1991 Snell et al 



19961. 



Ceccarelli et al. 



2000[lNisini et al.ll2002[ 



Maret et al. 



199i 



iarwit et al.lll99 



2002 



Boonman et al. 



Ceccarelli et al. 



20031. 



van der Tak et al 



20061 ) . Thus, water acts like a 'switch' that turns on whenever energy is deposited in molec- 
ular clouds and elucidates key episodes in the process of stellar birth, in particular when 
the system exchanges matter and energy with its environment. This includes basic stages 
like gravitational collapse, outflow injection, and stellar heating of disks and envelopes. Its 
unique ability to act as a natural filter for warm gas and probe cold gas in absorption make 
water a highly complementary diagnostic to the commonly used CO molecule. 

Studying water is also central to understandin g the fundamental proc esses of freeze- 



out, grain surface chemistry, and evaporation (e.g., iHoUenbach et al.l 120091 ). Water is the 



dominant ice cons tituent and can trap various molecules inside its matrix, including complex 



organic ones ( e.g.. lGibb et al.ll2004bl iBoogert et al.ll2008[ ). When and where water evaporates 
back into the gas is therefore relevant f or understanding the wealth o f organic molecules 
observed near protostars (see review by iHerbst fc van Dishoeckl |2009| ). Young stars also 
emit copious UV radiation and X-rays which affect the physical structu re as well as the 
chemistry, in particular that of hydrides like water (jStauber et al.l 120061 ). Moreover, the 
level of deuteration of water provides an important record of the temperature history of the 
cloud and the conditions during grain surface formation, and, in comparis on with cornetary 



data, of its evolution from interstellar cloud s to solar system o bjects (e.g.. IJacq et al.lll990 
Gensheimer et al.l Il996l . iHelmich et al.l Il996l iParise et al.l |2003[ ) . 



Finally, water plays aii actiy e role in the energy bala nce of dense gas (e.g.. lGoldsmith fc Langer 



19781 . iNeufeld fc KaufmanI Il993l . iDoty fc Neufeldl Il997l ). Because water has a large dipole 
moment, its emission lines can be efficient coolants of the gas, contributing significantly over 
the range of physical conditions appropriate in star-forming regions. The large dipole mo- 
ment can also lead to heating throug h absorption of infrared radiation followed by coUisional 
de-excitation (ICeccarelli et al.lll996l ). It is therefore important to study the gaseous water 
emission and absorption, and compare its cooling or heating efficiency quantitatively with 
that of other species. 



Interstellar water was detected more than 40 years ago through its 22 GHz maser emis- 
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sion towards Orion ( ICheung et al.lll969[ ). While this hne remains a useful beacon of star- 
formation activity, its special excitation and line formation conditions limit its usefulness as 
a quantitative physical and chemical tool. Because observations of thermally excited wa- 
ter lines from Earth are limited, most informati on on water has corn e from satellites. The 
Subm illimeter Wave Astronomy Satellite (SWAS. iMelnick et al.ll2000l ) and the Odin satellite 
(e.g., iHjalmarson et al.ll2003l . iBjerkeli et al.ll2009l ) observed only the 557 GHz ground-state 
line of ortho-H20 at spatial resolutions of ~ 3.3' x 4.5' and 126", respectively. The 557 
GHz line was detected in just the brightest objects by these missions and other lines of 
water (including those of para-H20) and related species could not be observed. The Short 
(SWS) and Long Wavelength Spectrometers (LWS) on the Infrared Space Observatory (ISO) 
covered a large number of pure rotational and vibration-rotation lines providing important 
insight in the water excitation and distribution, but with limited spe ctral and spatial res- 
olution and mapping capabilities (see review by Ivan Dishoeckl l2004j ) . The Spitzer Space 



Telescope ha s detected high 



from shocks (IMelnick et al. 



y excited pure rotationa. 



2008 



protoplanetary disks (ISalyk et al. 



H2O lines at mid- infrared wavelengths 
Watson et al.l 12007 ) and from the inner few AU reg ions of 



2008 



Carr fc Naiitall2008l . |Pontoppidan et al.ll2010[ ). Even 



higher-lying vibr ation-rotation H^O lines have been observed at near-infrared wavelengths 
from the ground (ISalyk et al.ll2008l ). However, these data do not provide any information on 
the cooler water reservoir where the bulk of the disk mass is. 

The 3.5m Herschel Space Observatory with, its suite of instruments ( Pilbratt et al. 2010 1^ 
is particularly well suited to address the distribution of cold and warm water in star- and 
planet-forming regions, building on the pioneering results from previous missions. Its wave- 
length coverage of 55-671 /im (0.45-5.4 THz; 15-180 cm~^ includes both low- and high- 
excitation lines of water, enabling a detailed analysis of its excitation and abundance struc- 
ture (Fig. [1]). Compared with the ISO-LWS beam (~80"), Herschel has a gain of a factor 
of ~ 8 in spatial resolution at similar wavelengths, and up to 10^ in spectral resolution. 
Compared to observatories with high spectral resolution heterodyne instruments, its diffrac- 
tion limited beam of 37" at 557 GHz is a factor of 3-5 smaller than that of SWAS or Odin, 
with a gain of > 10 in sensitivity because of the bigger dish a nd improved detector tech- 
nology. The Heterodyne Instrument for the Far- Infrared (HIFI; Ide Graauw et al.ll2010l ) has 
spectral resolving powers R = A/AA > 10^ up to 1900 GHz, allowing the kinematics of the 
water lines to be studied. The Photoconducting Array Camera and Spectrometer (PACS; 



^Herschel is an ESA space observatory with science instruments provided by European- led Principal 
Investigator consortia and with important participation from NASA 

^This paper follows the common usage of giving frequencies in GHz for lines observed with HIFI and 
wavelengths in fim for lines observed with PACS 
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Poglitsch et al.ll2010[ ) 5x5 pixel array provides instantaneous spectral mapping capabilities 
at i? = 1500 — 4000 of important (backbone) water lines in the 55-200 fim range. These 
enormous advances in sensitivity, spatial and spectral resolution, and wavelength coverage 
provide a unique opportunity to observe both cold and hot water in space, with no other 
space mission with similar capabilities being planned. 

The goal of the 'Water In Star-forming regions with Herschel' (WISH) Key Program 
(KP) is to use gas-phase water as a physical and chemical diagnostic and follow its abundance 
throughout the different phases of star and planet formation. A comprehensive set of water 
observations is carried out with HIFI and PACS toward a large sample of young stellar 
objects (YSOs), covering a wide range of masses and luminosities - from the lowest to the 
highest mass protostars- and a large range of evolutionary stages - from the first stages 
represented by the pre-stellar cores to the late stages represented by the pre-main sequence 
stars surrounded only by disks. Lines of H2O, Hl^O, Hg^O and the chemically related species 
O, OH, and HaO^ are targeteclfl. In addition, a number of hydrides which are diagnostic 
of the presence of X-rays and UV radiation are observed. Selected high-frequency lines of 
CO, ^^CO, and C^^O as well as dust continuum maps are obtained to constrain the physical 
structure of the sources independent of the chemical effects. Together with the atomic O 
and C"*" lines, these data also constrain the contributions from the major coolants. The 
Herschel data are complemented by ground-based maps of longer wavelength continuum 
emission and lines of HDO, CO, C and other molecules to ensure a self-consistent data set 
for analysis. In terms of water, the WI SH program is comple mentary to the 'CHEmical 



Survey of Star-forming regions' f CH ESS: ICeccare. 



of Extraordinary Sources' (HEXOS; iBergin et al 



li et al.ll2010[ ) and 'Herschel observations 



2010bl ) HIFI KPs which survey the entire 



spectral range (including many lines of water) but only for a limited number of sources 
and with smaller on-source integration times per frequency setting. The 'Dust, Ice and 
Gas in Time' (DIGIT) (PI N.J. Evans) and 'Herschel Orion Protostar Survey' (HOPS) 
(PI T. Megeath) KPs complement WISH by ca rrying out full PACS spe ctral scans for a 



larger sample of low-mass embedded YSOs (e.g.. Ivan Kempen et al.ll2010a] ). The 'PRobing 
Interstellar Molecules with A bsorption line Stu dies' (PRISMAS) targets the water chemistry 
in the diffuse interstellar gas (iGerin et al.ll2010l ) whereas 'W ater and related chem istry in the 
Solar System' (HssO) observes water in planets and comets (IHartogh et al.ll2010l ). Together, 
these Herschel data on water and related hydride lines will provide a legacy for decades to 



come. 



In the following sections, background information on the water chemistry and excitation 



^Both H2O and 'water' are used to denote the main Hj^O isotopologue throughout this paper 
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relevant for interpreting the data is provided (§2). The observational details and organization 
of the WISH program are subsequently described in §3, with the specific goals and first results 
of the various subprograms presented in §4. A discussion of the results across the various 
evolutionary stages and luminosities is contained in §5 together with implications for the 
water chemistry, with conclusions in §6. Detailed information can be found at the WISH 
KP website http : / /www . strw . leidenuniv . nl/WISH/ , This website includes links to model 
results, modeling tools, complementary data and outreach material. These data and analysis 
tools will be useful not only for Herschel but also for planning observations with the Atacama 
Large Millimeter /sub millimeter Array (ALMA) and future far-infrared missions. 



2. H2O chemistry and excitation 

Star formation takes place in cold dense cores with temperatures around 10 K and H2 
densities of at least 10^ cm~'^. Once collapse starts and a protostellar object has formed, 
its central luminosity heats the surrounding envelope to temperatures well above 100 K. 
Moreover, bipolar jets and winds interact with the envelope and cloud, creating shocks 
in which the temperature is raised to several thousand K. Cloud core rotation leads to a 
circumstellar disk around the young star with densities in the midplane of > 10^° cm^'^. 
With time, the envelope is dispersed by the action of the outflow, leaving the pre-main 
sequence star with a disk only. Once accretion onto the disk stops and the disk becomes 
less turbulent, the ~0.1 /im grains from the collapsing cloud coagulate to larger particles 
and settle to the midplane, eventually leading to planetesimals and (proto)planets. Thus, 
star-forming regions contain gas with a large range of temperatures and densities, from 10 to 
2000 K and 10^ to > 10^*^ cm~^, and with gas/dust ratios that may vary from the canonical 
value of 100 to much larger or smaller values. The water chemistry responds to these different 
conditions. 



2.1. Chemistry 



Figure [2] contains an overview of the physical structure of a typical protostellar en- 
velope, with the main chemical regimes indicated. In standard gas-phase chemistry, H2O 
forms through ion-molecule reactions sta r ting with O + Ht and O'*' + H2, both leading 
to 0H+ (e.g., iHerbst &: Klempererl 1 19731 . iDalgarno fc Blackl flQTel ) (see Fig. [3]). A series 
of rapid hydrogen abstraction reactions with H2 leads to HaO"*", which can dissociatively 
recombine to form H^O and OH with branching ratios of ~33% and 66%, respectively 
(IVejby-Christensen et al.l 119971 ). H2O is destroyed by ultraviolet photons and by reactions 
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Fig. 1. — Energy levels of ortlio- and para-H20, with HIFI transitions (in GHz, red) and 
PACS transitions (in /xm, blue) observed in WISH indicated. 




Fig. 2. — Temperature (blue) and density (black) structure of a high-mass protostellar 
envelope. For illustrative purposes, a typical abundance profile is included (red dash-dotted 
line) on an arbitrary scale, with the relevant chemical processes in each regime indicated 
at the top of the figure. The figure is ba sed on the model of AFGL 2591 (L = 6 x 10^ 
Lq) presented by Ivan der Tak et al.l ( 119991 ). For distances of ~1 kpc, the typical Herschel 
beam of 20" corresponds to 20000 AU, i.e., it encompasses the entire envelope. For low- mass 
protostars, the ice evaporation zone moves inward to radii typically less than 100 AU. 
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Fig. 3. — Summary of the main gas-phase and sohd-state chemical reactions leading to the 
formation and destruction of H2O. Three different chemical regimes can be distinguished: (i) 
ion-molecule chemistry, which dominates gas-phase chemistry at low T; (ii) high-temperature 
neutral-neutral chemistry; and (iii) so lid state chernistry. The latter chemical scheme is 
based on the latest laboratory data by lloppolo et al.l (l2010l ). s-X denotes species X on the 
ice surfaces. 
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with C 



Hg and other ions such as HCO 



All key reactions in this netw ork have been mea- 



sured or calculated at low temperatures (e.g.. iKlippenstein et al.ll2010l ). The steady-state 
H2O abundance resulting from pure gas-phase chemistry is typically a few ti mes 10 ~^ with 
respe ct to H2 at low T and scales with the cosmic-ray ionization rate ( (e.g. 

mm- 



Woodall et al. 



In the last decades, this traditional view of the oxygen chemistry has been radically 
changed, triggered by the very low abundances of O2 in cold clouds found by SWAS (e.g., 
Bergin et al.l 120001 ). It is now realized that in cold and dense regions, H2O is formed much 
more efficiently on the grains through a series of reactions involving O and H accreted 
from the gas, as predi cted more than 50 years ago by van de Hulst and quantified by 
Tielens &: HagenI (jl982l ) . Surface s cience laboratory experiments are only now starting to 
test these reactions experimentally (lloppolo et al.ll2008l . Iloppolo et al.ll2010l . iMiyauchi et al. 



2008, Duheu et al. 



2010, 



stood (e.g.. lCuppen et al. 



Romanzin et al. 201ol ). While the details are not yet fully under- 
2OI0I ). there is ample observational evidence from 3 /^m absorption 



spectroscopy of the solid H2O band toward YSOs that the H2O ice abundances can be as high 



as ~ 10 ^, much higher than results from freeze-out of gas-phase H2O (e.g., IWhittet et al. 



19881 . IPontoppidan et al.ll2004[ ). 



Water ice f ormation start s at a threshold extinction of Ay ~ 5 mag, i.e., well before the 



cloud collapses (IWhittetll2003l ). A small fraction of the H2O ice formed by surface reactions 
can be desorbed back into the gas phase by a variety of non-thermal pr ocesses including 



cosmic-ray induced desorption and photodesorption (see recent overview by lHoUenbach et al. 



2OO9I ) . The efficiency of the latter p rocess has been quantified by laborator y experiments to 
be ~ 10-3 per incident UV photon JWestlev et al.lll995l . [6berg et al.ll2009h . 



Near protostars, the grain ter nperature rises abo ve ~ 100 K and all the H2O ice thermally 
desorbs on very short timescales (IFraser et al.ll200ll ). leading to initial gas-phase abundances 
of II2O as high as the original ice abundances of ~ 10"^. Other molecules trapped in the 
water ice matrix evaporate at the same time in these so-called 'hot core' or 'hot corino' (in 
the case of low-mass YSOs) regions. Thus, constraining the H2O evaporation zone is also 
critical for the interpretation of complex organic molecules in hot cores. 

At even higher temperatures, above 230 K, the gas -phase reaction O + Ho — )• OH + 
H, w hich is endoergic by ~2000 K, becomes significant (lElitzur fc Watson! 1 19781 . ICharnley 
19971 ). OH subsequently reacts wi th H^ to form H3O, a reaction which is exothermic but 
has an energy barrier of ~2100 K (lAtkinson et al.ll2004j ). This route drives all the available 
gas-phase oxygen into H2O leading to an abundance of ~ 3 x 10""^ at high temperatures, 
unless the H/H2 ratio in the gas is so high that the back-reactions become important. Such 
hot gas can be found in the hot cores close to the protostars themselves and in the more 
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extended shocked gas associated with the outflows. Given the range of temperatures and 
densities in protostellar (Fig. |2]) and protoplanetary environments, all of the above chemical 
processes likely play a role, and the Herschel data are needed to determine their relative 
importance. 

H2O has a unique deuteration pattern compared with other molecules. If grain surface 
formation dominates, water is deuterated in the ice at a level that may be orders of mag- 



nitude lower than that of other s 



Decies which participate more actively in the dense cold 



gas chemistry (IRoberts et al.ll2003l ). Indeed, t he deuteration fraction observed in high mas s 



[e.g. 



Gensheimer et al.lll996l iHelmich et al. 1996) 



hot cores is typically HD0/H20< 10 

although higher val ues around 10~^ have recently been found for Orion (iPersson et al 



12007 



Bergin et al.ll2010al ) and for a low- mass hot core (jParise et al.l l2005bl ) . All of these values 



are still much lower th an, for examp 



e, DCN/HCN, HDCO/H2CO fe.g.. iLoinard et aimOlh 



or CH2DOH/CH3OH (iParise et al.ll2004f). I n the sohd phase, the HDO/H2O upper limits 
are less than 2 x 10"'^ ( Dartois et al. 2003 . Parise et al. 20031 ). An important question is 



how similar the protostellar HDO/H2O ratios are to those observed in comets and in water 
in our oceans on Earth (about 2 x 10~^), since this has major implications for the del ivery 
mechanisms of water to our planet (IBockelee-Morvan et al.lll998l . [Raymond et al.ll2004l ). No 
specific HDO lines are targeted with H erschel within the W ISH program, since several HDO 
lines can be observed from the ground (jParise et al.ll2005a[ ). HIFI will provide the necessary 
data on H2O itself for comparison with those HDO data. 



2.2. Excitation 

H2O is an asymmetric rotor with a highly irregular set of energy levels characterized by 
quantum numbers JkaKc ■ Because of the nuclear spin statistics of the two hydrogen atoms, 
the H2O energy levels are grouped into ortho {Ka + Kc =odd) and para {Ka + Kc =even) 
ladders across which no transitions can occur except through chemical reactions which ex- 
change a H nucleus (Fig. [T]). The H2O energy levels within each ladder are populated by 
a combination of collisional and radiative processes. The most important collision partners 
are ortho- and para-H2, with electrons and H only significant in specific regions and He 
contributing at a low level. Accurate collisional rate coefficients are thus crucial to interpret 
the Herschel data. The bulk of these rates are obtained from theoretical quantum chemistry 
calculations which involve two steps. First, a mult i- dimensional potential energy surface 
involving the colliders is computed. Second, the dynamics on this surface are investigated 
with molecular scattering calculations at a ra nge of collision eii ergies. Early studies used 



simplifications such as replacing H2 with He (IGreen et al.lll993l ). including only a limited 
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number of degrees of freedom of the co lliding system, or using approximations in the scat- 
tering calculations (jPhillips et al.lll996l ). Recent calculations consider a 5D potential surface 
and compute collisions with both para-H^^ and ortho-Ho separately up to h igh temperatures 
with the close coupling method (jPubernet et al.ll2009l . iDariiel et al.l 120101 ) . Thi s 5D poten- 
tial is based on the full 9D potential surface computed by IValiron et al.l ( l2008l ). The cross 
sections for collisions of H2O with 0—B.2 are found to be significantly larger than those with 
|}-H2 for low J, so that explicit treatment of the o/p H2 ratio is important. 

Direct comparison of absolute values of theoretical state-to-state cross sections with 
experimental data is not possible, but the accuracy of the 9D potential surface has been 
confirmed implicitly through comp arison with other data sets, including differential scatter- 
ing experiments of H2O with H2 (lYang et al.ll2010l ). An indication of the accu racy of the 



calculations can also be obtained from recent pressure broadening experiments (jPick et al. 



20101). The high temperature (> 80 K) data agree with theory within 30% or better, but the 
low temperature experimental rate coefficients are an order of magnitude smaller. This dis- 
crepancy has recently been understood by realizing that the classical impact approximation 
used in the analysis of the pressure broadening data does not hold at very l ow temperatures 



( Wiesenfeld fc Faurd l2010l ). The computed rate coefficients at low T by iDubernet et al. 
J2006h should therefore be valid. 



Because of the large rotation constants of water, the energy spacing between the lower 
levels is much larger than that of a heavy molecule like CO. Thus, the water transitions 
couple much more efficien tly with far-infrared r adiation from warm dust, which can pump 
higher energy levels (e.g., iTakahashi et al.l 119831 ). This far- infrared continuum also plays a 
role in the line formation (e.g., whether the line occurs in absorption or emission) and thus 
provides constraints on source geometry. Moreover, the continuum can become optically 
thick at the highest freque ncies, providing an effective 'screen' for looking to differ ent depths 
in the YSO environment (jPoelman fc van der TakI 120071 . Ivan Kempen et al.ll2008l ). 



The high transition frequencies combined with the large dipole moment of H2O (1.85 
Debye) also lead to high line opacities. For example, for the o— H2O 557 GHz ground state 
transition, the line center optical depth is tq = 2.0 x 10^^^N£{o — H2O) A\^, where Ni is the 
column density in the lower level in cm~^, which is to first approximatio n equal to the tota l 
o— H2O column density, and AV is the FHWM of the line in km s~^ ( jPlume et al.ll2004j ). 



Thus, even for H2O abundances as low as 10 in typical clouds with A^(H2)> 10 cm , 
the water lines are optically thick at line center. 

Because of the high optically depths, it is difficult to extract reliable water abundances 
from lines of the main isotopologue. However, because the excitation is subthermal (critical 
densities are typically 10^—10^ cm~^), every photon that is created will scatter and eventually 
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escape the cloud. In this so-called effectively thin limit, the water column density scales 
linearly with integrated intensity of a ground-state transitio n and is inversely proportional 



to de nsity and the value of the coUisional rate coefficient (ISnell et al.l l2000l ISchulz et al. 
199lh . 



Many of the H2O transitions exhibit population inversion, either very weakly or more 
strongly, suc h as the famous 61 r — 823 maser at 22 GHz widely observed in star- forming 
regions (e.g., ISzymczak et al.l 120051 ). Thus, the H2O molecule forms a formidable challenge 
for radiative transfer codes, and the HIFI and PACS spectra provide a unique reference data 
set of thermal lines against which to test the basic assumptions in the maser models. 



2.3. Modeling tools 



The WISH team has developed a variety of tools important for the Herschel data anal- 
ysis, several of which are publicly available to the community. This includes a molecular 
line database LAMDA0 for excitation calculations of H2O, O H, CO and other m olecules 
using various sets of coUisional rate coefficients discussed above ( ISchoier et al.ll2005[ ). Simple 
ID escape probability non-LTE radiative transfer programs such as R ADEX can be run 



eithe r on-line or downloaded to run at the researcher's own institute (Ivan der Tak et al. 



20071 ). This and the more sophisticated Monte-Carlo ID radiative transfer code RATRAN 
are publicly availabl^ and RADEX is also included in the CASSIS set of analysis pro- 
gramj^. The RATRAN c ode has been extensively tested in a code comparison campaign 

3 ' 



van Zadelhoff et al.ll2002r) . Another ID non-LTE code used by the WISH team is MOLLIE 



( iKeto et al.ll2004j ). A 2D version of RATRAN as well as a ne w flexible 3D line radiative 



trans f er code called LIME are avai lable on a collaborative basis (IHogerheijde fc van der Tak 
2OOOI. iBrinch &: Hqgerheiidd I2OIOI ). A 2D escape probability code has been written by 



been developed by 



Poelman fc Spaans (l2005l) wherea s a fast 3D code using a local source approximation has 



Bruderer et al.l ( l2010al ) 



^http : //www . strw . leidenuniv . nl/ ^moldata 

" http : / / www . sron . rug . nl/~vdtak/ ratrein/f reunes . html 



http: //www. cesr.fr/~walters/web_cassis/ 
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Table 1. WISH Source List 



Coordinates Properties 



RA Dec Vlsr ^bol d Outflow'' RD'' Ref. Notes 

(hms) (°"') (kms-l) (Lq) (pc) (" ") 



Prestellar cores 



L 1544 05 04 17.2 +25 10 42.8 +7.3 • • • 140 

B68 17 22 38.2 -23 49 54.0 +3.4 ••• 125 



Low-mass YSOs 





UO 


9"^ 


oc5.y 


1. on 




A 


1 c q 


1 1 fi 
11. D 


9 Pin 


T\ ( ^^ -U9J^ 

-D — 10, -\-^0.<J) 




1 9 


iN^jrO looo irtiVo Z 




OQ 
ZO 


uo.o 


1 Q1 

-hoi 


1 A 


'i'7 1 
o( .1 




on v 


Zoo 


13 — yz, -\-zo) 

■R ( -i.7(\ 1 K\ 


+ 


4,2, d 




03 


29 


03.8 


-roi 


16 


04.0 


+8.5 


50 


235 


JJX I T^ZUj — ZU J 

B2 (+20, -50) 


+ 


1,0,0 


NGC1333 IRAS4A 


03 


29 


10.5 


+31 


13 


30.9 


+7.2 


7.7 


235 


B (-6, -19) 
R (+13, +25) 




1,2,3 


NGC1333 IRAS4B 


03 


29 


12.0 


+31 


13 


08.1 


+7.4 


7.7 


235 






1,2,3 


L1527 


04 


39 


53.9 


+26 


03 


09.8 


+5.9 


2 


140 


B (+40, +10) 




1,2,6 


BHR71 


12 


01 


36.3 


-65 


08 


53.0 


-4.4 


10 


200 


B (+40, -100) 
R (-40, +140) 




6,7,8 


L 483 MM 


18 


17 


29.9 


-04 


39 


39.5 


+5.2 


9 


200 


B (-40, 0) 




1,2,6 


Scr SMM 1 


18 


29 


49.8 


+01 


15 


20.5 


+8.5 


30 


250' 


Bl (-15, +30) 


+ 


9,10,11 


Ser SMM 4 


18 


29 


56.6 


+01 


13 


15.1 


+8.0 


5.0 


250' 


R (+30, -60) 




9,12,11 


Ser SMM 3 


18 


29 


59.2 


+01 


14 


00.3 


+7.6 


5.9 


250' 






9,12,11 


L723 MM 


19 


17 


53.7 


+19 


12 


20.0 


+11.2 


3 


300 






6,13 


B335 


19 


37 


00.9 


+07 


34 


09.6 


+8.4 


3 


250 


B (+30, 0) 




1,2,6 


L1157 


20 


39 


06.3 


+68 


02 


15.8 


+2.6 


6 


325 


B2 (+35, -95) 
R (-30, +140) 




1,2,6 


L1489 


04 


04 


43.0 


+26 


18 


57.0 


+7.2 


3.7 


140 




+ 


14,15,6 


L 1551 IRS 5 


04 


31 


34.1 


+18 


08 


05.0 


+7.2 


28 


140 


B (-255, -255) 
R (+150, +20) 




16,12,6 


TMRl™ 


04 


39 


13.7 


+25 


53 


21.0 


+6.3 


3.7 


140 






4,12,6 


TMC lA™ 


04 


39 


34.9 


+25 


41 


45.0 


+6.6 


2.2 


140 






12,6 


TMCl 


04 


41 


12.4 


+25 


46 


36.0 


+5.2 


0.7 


140 






12,6 


HH46 


08 


25 


43.9 


-51 


00 


36.0 


+5.2 


12 


450 


B (-10, 0) 
R (-40, -20) 




17,18 


CedllO IRS4 


11 


06 


47.0 


-77 


22 


32.4 


+3.5 


1 


125 






19,1,20 


IRAS 12496/DK Cha 


12 


53 


17.2 


-77 


07 


10.6 


+2.3 


50 


200 






1,20 


IRAS 15398/B228 


15 


43 


01.3 


-34 


09 


15.0 


+5.1 


1 


130 






1,20 


GSS 30 IRSl 


16 


26 


21.4 


-24 


23 


04.0 


+2.8 


25 


125 






1,4,21 


Ellas 29 


16 


27 


09.4 


-24 


37 


19.6 


+5.0 


36 


125 






1,22,21 


Oph IRS 63 


16 


31 


35.6 


-24 


01 


29.6 


+2.8 


1.6 


125 






1,22,21 
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Table 1 — Continued 





Coordinates 




Properties 








RA 


Dec 


Vlsr 


Lhoi 


d 


Outflow'^ 


RD'' Ref. Notes 




(h m s) 




(kms-i) 


(Le) 


(pc) 


(" ") 




RNO 91 


16 34 29.3 


-15 47 01.4 


+5.0 


11 


125 




23,1,21 


R CrA IRS 5 


19 01 48.0 


-36 57 21.6 


+5.7 


3 


170 




24,1,20 


HH 100 


19 01 49.1 


-36 58 16.0 


+5.6 


14 


125 




25,1,20 


Intermediate-mass YSOs 


AFGL490 


03 27 38.4 


+58 47 08.0 


-13.5 


2000 


1000 




26,27,28 


L1641 S3MMS1 


05 39 55.9 


-07 30 28.0 


+5.3 


70 


465 




29,30 


NGC 2071 


05 47 04.4 


+00 21 49.0 


+9.6 


520 


422 


PO (+32, +63) 


31,32,30 














P5 (-128, -97) 




Vela IRS 17 


08 46 34.7 


-43 54 30.5 


+3.9 


715 


700 




33,34,35 


Vela IRS 19 


08 48 48.5 


-45 32 29.0 


+12.2 


776 


700 




33,35 


NGC 7129 FIRS 2 


21 43 01.7 


+66 03 23.6 


-9.8 


430 


1250 


B (+60, +60) 


+ 36,37,38 














R (+60, -60) 




Outflow only 


HH211-mm 


03 43 56.8 


+32 00 50 


+9.2 


4 


250 


C (0, 0) 
















B2 (+37, -15) 




IRAS 04166 


04 19 42.6 


+27 13 38 


+6.7 


0.4 


140 


B (+20, +35) 
















R (-20, -35) 




VLA-1 HHl/2 


05 36 22.8 


-06 46 07 


+8 


50 


450 


B (+60, -80) 




HH212 MMl 


05 43 51.4 


-01 02 53 


+1.6 


15 


460 


C (0, 0) 
















B (-15, -35) 




HH25 MMS 


05 46 07.3 


-00 13 30 


+10 


6 


400 


C (0, 0) 
















SiO (+36, -57) 




HHlll VLAl 


05 51 46.3 


+02 48 30 


+8.5 


25 


460 


C (0, 0) 
















Bl (-170, +21) 




HH54B 


12 55 50.3 


-76 56 23 


+2.4 


1 


180 


C (0, 0) 




VLA 1623 


16 26 26.4 


-24 24 30 


+3.5 


1 


125 


Bl (+30, -20) 
















R2 (-65, +25) 




IRAS 16293 


16 32 22.8 


-24 28 36 


+4.5 


14 


125 


B (+75, -60) 
















R (+75, +45) 




Ser S 68 N 


18 29 47.5 


+01 16 51 


+8.5 


4 


260 


B (-12, +24) 
















C (0, 0) 




Cep E MM 


23 03 13.1 


+61 42 26 


-11.0 


75 


730 


B (-10, -20) 










High-mass 


YSOs 








GII.II-O.I2-NH3-PI 


18 10 33.9 


-19 21 48 


+30.4 




3600 




39, 40 « 


G11.11-0.12-SCUBA-P1 


18 10 28.4 


-19 22 29 


+29.2 




3600 




40 


G28.34+0.06-NH3-P3 


18 42 46.4 


-04 04 12 


+80.2 




4800 




40 



-17- 



Table 1 — Continued 



Coordinates Properties 





RA 


Dec 


^LSR 




d 


Outflow^ 


RD^ Ref. 


Notes 




(h m s) 


(O , 


(kms-l) 




(pc) 


(" ") 






G28.34+0.06-SCUBA-P2 


18 42 52.4 


-03 59 54 


+78.5 




4800 




40 




IRAS05358+3543 


05 39 13.1 


+35 45 50 


-17.6 


6.3 X 10=^ 


1800 




41 


f 


IRAS16272-4837 


16 30 58.7 


-48 43 55 


-46.2 


2.4 X lO'' 


3400 




42 




NGC6334-I 


17 20 53.3 


-35 47 00 


-4.5 


1.7 X 10"* 


1700 




+ 43 




W43-MM1 


18 47 47.0 


-01 54 28 


+98.8 


2.3 X 10* 


5500 




44 




DR21(OH) 


20 39 00.8 


+42 22 48 


-4.5 


1.7 X lO"* 


1700 




45, 46 




W3-IRS5 


02 25 40.6 


+62 05 51 


-38.4 


1.7 X 10"' 


2200 




+ 47, 48, 49 


g 


IRAS18089-1732 


18 11 51.5 


-17 31 29 


+33.8 


3.2 X 10" 


3600 




50 




W33A 


18 14 39.5 


-17 52 00 


+37.5 


1.0 X 10"* 


4000 




51 




IRAS18151-1208 


18 17 58 


-12 07 27 


+32.0 


2.0 X 10* 


2900 




52 




AFGL2591 


20 29 24.9 


+40 11 19.5 


-5.5 


5.8 X 10* 


1700 




+ 53 




G327-0.6 


15 53 08.8 


-54 37 01 


-45.0 


1.0 X 10* 


3000 




54, 55 


h 


NGC6334-I(N) 


17 20 55.2 


-35 45 04 


-7.7 


1.1 X 10^ 


1700 




+ 56 




G29.96-0.02 


18 46 03.8 


-02 39 22 


+98.7 


1.2 X 10^ 


7400 




57, 58 




G31.41+0.31 


18 47 34.3 


-01 12 46 


+98.8 


1.8 X lO"' 


7900 




59, 60, 61 




IRAS20126+4104 


20 14 25.1 


+41 13 32 


-3.8 


1.0 X 10* 


1700 




50 




G5.89-0.39 


18 00 30.4 


-24 04 02 


+10.0 


2.5 X 10* 


2000 




62, 63 


1 


GlO.47+0.03 


18 08 38.2 


-19 51 50 


+67.0 


1.1 X 10* 


5800 




64, 65, 57 




G34.26+0.15 


18 53 18.6 


+01 14 58 


+57.2 


2.8 X 10* 


3300 




66, 65 




W51N-C1 


19 23 43.8 


+14 30 26 


+59.5 


1 - 10 X lO'^ 


5500 




62 




NGC7538-IRS1 


23 13 45.3 


+61 28 10 


-57.4 


2.0 X 10* 


2800 




+ 67, 68 




Radiation diagnostics 


S140 


22 19 18.2 


+63 18 46.9 


-7.1 


10* 


910 








Disks 


DM Tau 


04;33;49.7 


+18:10:10 


6.1 


0.25 


140 






j 


LKCalS 


04:39:17.8 


+22:21:04 


6.0 


0.74 


140 








MWC480 


04:58:46.3 


+29:50:37 


5.0 


11.5 


140 








TW Hya 


11:01:51.9 


-34:42:17 


3.0 


0.25 


56 








BP Tau 


04:19:15.8 


+29:06:27 


6.0 


0.9 


140 






k 


GG Tau 


04:32:30.3 


+17:31:41 


7.0 


0.8, 0.7 


140 








GM Aur 


04:55:10.2 


+30:21:58 


5.0 


0.7 


140 








MWC758 


05:30:27.5 


+25:19:57 


5.5 


21. 


200 








T Cha 


11:57:13.5 


-79:21:32 


5.0 


0.4 


66 








IM Lup 


15:56:09.2 


-37:56:06 


4.5 


1.7 


150 








AS209 


16:49:15.3 


-14:22:07 


6.6 


0.7 


160 








HD163296 


17:56:21.3 


-21:57:20 


6.0 


35. 


122 









Also part of the Outflow sub-program, and will be observed at the given offset coordinates 
Also part of the Radiation Diagnostics sub-program 
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Low-mass class YSOs 

Low-mass class I YSOs 

High-mass prestellar cores 
f mIR-quiet HMPOs 
s mIR-bright HMPOs 
^ Hot Molecular Cores 
' UC HH Regions 
j Gas-rich Disks, deep sample 
^ Gas-rich Disks, shallow sample 

' Using VLBA observations of a star thought to be associated with the Serpens cluster, iDzib et al. I ll20ld) derive a distance of 415 



pc 



™ Coordinates used in WISH; the mo re accurate coordinate s obtained from SMA observations are 04 39 13.9, -1-25 53 20.6 (TMRl) 
and 04 39 35.2, +25 41 44.4 (TMCIA) l|j0reensen et al.ll2009h 



References. — (1 ) lEvans et al] i200^): (2) Ijprgensen et al.l ll2007h : (3) Iffirota et al.l ||2008|'): (4) lj0rgensen et al.l ll2009h: (5) 



iBachiller et al.l lll998l'); (6) Andre ct al. (2000); (7) Chen ct al.' ('2008''); fS) ISeidenstick er fc Schmidt-Kalcr' ('l989^: (91 iKristensen et al. 



I); (10) Ivan Kempen et al.^_ (2009c'); (1 1) Pgntoppid a ,n ct a lj (l2004h ; (12 
14) iHogerheiide et al.l (1199?!); (15) iBrinch et al.[ i2007f); 



iNoriega-Crespo et al. ( 20041 ); (19) 



IgSeJ); (28) I Snell etaU |[l98^ 
1 1990-) 



11998) 



Knude , 



..Froebrich[(l2005l ); (20) 
(23) | Chen et al.l ||200£ ) ; (24) iNisini et al.l i 2005^; (25) Iwilking et al. 



j2( 



(29) LStanke et al.' (2000); (30) 'Wilson ct al. 



LI16UBU 

Hod mm . 

'aiT(ll992[) ; 



Butner et al, 

(21)lde Geus et 



j0rgensen et al. 



(26) ICuUixson 
2005'): (31) 'Johnstone et al. 



2007); (18) 



Hoi^erheiide et al.. (1999.') ; (13) 

(Il99ll); (17) l yelusamv et al . . _ 

al] lll990h: ( 22) lLommen et al.ril2008h: 
et al.ri983|) ; (27) IMozurkewich et al ' 



(2001); (32) Butner et al. 



(33) Lis eau et al]lll992l): (34) [Giannini ct al. (2005); (35) Slawson & Reed (1988); (36) Fucntc ct al. (2005a); (37) Eiro a et al, 
(38) She vchenko fc Yakubovl Il989l): (39) ICarev et al.l ||1998D: (40) IPillai et al.l j20l36bi ): (41) iBeuther et al.l (l2 002a): (42) 



Gara^_ot_aL 



(43) 



Beuther et al 



(I2OO8I' 



\k (44) iMotte et al. 



(52) [ Beuther et al.l (120021: ) 



Beuther et al 



jvarT'der Tak etal 
(58) iGibb et al.l 



2001 



(12004: 



_12( 

(45) Csengeri et al. in pre ss; (4 6)lMotte et al.l (t2007.): (47 ) 



Rodon et all (|200^ ): (48)lHachisuka et al.l (I2OO6I): (49 ) lLadd et al1 ( 1993l): (50)ISridharan et al.l ( l20o3): (51)lvan der Tak et alT ( 2000d) : 



(54 ') IWvrowski et al.l ll2006h : (55) iBronfman et al.l ( 1199^) ^ : (56) ISandell 
(5 9)lBeltran et al.l (20051 ): (60) ICesaroni et al.l (Il998l ): ( 61)lMueller et al. 



(57) ,_ . r- -1 

(62)'Somns ct al.'('2004'): (63)l Hunter et a l.'('2008): (64)'01mi et al.l (Il996l) : (65) "Wyrowski (prr 
(67) ,Sandell Sievers, (.2004, ): (68) lDe Buizer fc Minier, 12001 



iv. comm 



,); (66) lMookeri ea et al. 
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Recognizing the importance and complexity of the radiative transfer problem for H2O, 
the HIFI consortium organized a benchmark workshop dedicated to an accurate comparison 
of radiative transfer codes in 2004. The model tests and results are available at a Web 
pag^ so that future researchers can test new codes. The reliability of existing codes has 
significantly improved thanks to these efforts. 

Predictions of water emission lines have been made for grids of ID spherically symmetric 
envelope models over a large range of lumi nosities, envelope masses and other YSO parame 
ters, for a range of H2O abundances (e.g.. iPoelman fc van der Takll2007l . Ivan Kempen et al. 
20081 ). Many of the line profiles show deep absorptions and horn-like shapes due to the high 
line center optical depths discussed above, with emission only escaping in the line wings. 

Chemical mod els of envelopes a r id ho t cores have been developed for low- and high- 
mass sources (e.g.. 



Dotv et al. 


2002; 


2004, 


Viti et al. 


2001) 



cal conditions are kept static with time. Models in which the physical structure of the 
so urce changes with time a s the cloud collap s es an d the luminosity evolve s have been made 



Rodeers & Charnlev 


( 


200c 


), 


Lee et al. 


in ID and by 


Visser et al. ( 


2009) in 2D. 



Grids of ID C— type sh ock models by Kaufman fc Neufeldl (119961 ) and bow-shock models by 
Gustafsson et al.l ( 120101 ) are also available. The effect s of high-energy irra.diation have been 



studied in spherical symmetry in a series of papers (IStauber et al 



have been extended to multidimensional geometries (iBruderer et al 



2004; 



2005 



2009b 



a 



2010af ). 



20061) and 



2.4. Modeling approach 

Spherically symmetric models of protostellar envelopes such as illustrated in Fig. |2] are 
constructed for all sources in the WISH sample by assuming a power-law density profile and 
calculating t he dust temperature w ith radius usin g a continuum radiative trans fer code such 
as DUSTY Jivezic fc Elitzuil Il997h or RADMC jDuUemond fc DominikI I2004} ). taking the 
central luminosity as input. The power- law exponent, spatial extent and envelope dust mass 
are determined from minimization to the far- infrared spectral ener gy distribution and 
submillimeter continuum maps (see procedure by |j0rgensen et al.l 120021 ). The gas tempera- 
ture is taken to be equal to the dust temperature, wh ich is a good approxirn ation at these 



high densities where gas-dust coupling is significant (iDoty fc Neufeldl 1 19971 ). and the gas 
mass density is obtained through multiplication of the dust density by a factor of 100. These 
models are termed 'passively heated' envelope models, to distinguish them from 'active' 



^ttp : //www ■ sron . rug . nl/ ~vdtak/H20| 
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shock and UV photon heating mechanisms. 



At typical distances of star-forming regions, the warm and cold regions - and thus the 
H2O chemistry zones- will be largely spatially unresolved in the Herschel beams. However, 
the abundance variations ('where is the water spatially along the line of sight?') can be recon- 
structed through multi-line, single position observations coupled with physical models of the 
sources and radiative transfer analyses. This 'backward modeling' technique for retrieval of 



on a variety of molecules such as CO (, 


J0re:ensen et al 


200 


2), CH3C 


)H fe.e.. van der Tak et al. 


2000a. Maret et al. 2005 


Kristensen et al. 


2010a) and HoCO (e.g. 


van Dishoec 


: et al. 


I995I, 


CeccareUi et al. 2000bl. ! 


3choier et al. 


2002 


■ 


Vlaret et al. 


2005. J0re:ensen et al. 


2005b 


) and 


confirmed bv interferometer data (J0re;ensen 


2004. e 


khoier et al. 


2004). Asvmmetric rotors 



like H2O with many lines from different energies close in frequency are particularly well 
suited for such analyses. 



Table 2. WISH Line List 



Species Transition v A E^/k^B A Ref. 

(GHz) (Mm) (K) (s-i) 



HIFI lines 



H20 


lio-loi 


556.936 


538.3 


61.0 


3.5( 


-3) 


1 




2i2 — Iqi 


1669.905 


179.5 


114.4 


5.6( 


-2) 


1 






1113.343 


269.3 


53.4 


1.8( 


-2) 


1 




2o2 — 111 


987.927 


303.5 


100.8 


5.8{ 


-3) 


1 






752.033 


398.6 


136.9 


7.1{ 


-31 


1 




3i2-3ot 


1097.365 


273.2 


249.4 


1.6( 


-21 


1 




312-221 


1153.127 


260.0 


249.4 


2.6{ 


-3) 


1 


J.J.2 


lio-loi 


547.676 


547.4 


60.5 


3.3( 


-3) 


1 




111 — Ooo 


1101.698 


272.1 


52.9 


1.8( 


-2) 


1 






994.675 


301.4 


100.6 


6.0( 


-3) 


1 






1095.627 


273.6 


248.7 


1.6( 


-2) 


1 


J.J.2 vy 




552.021 


543.1 


60.7 


3.4( 


-3) 


1 




lii~Ou() 


1107.167 


270.8 


53.1 


1.8{ 


-2) 


1 




J =1/2,3/2-1/2,1/2 


1834.747 


163.4 


269.8 


6.4( 


-2) 


1 




f2, J =1/2,3/2-1/2,1/2 


1837.817 


163.1 


270.1 


6.4( 


-2) 


1 


OH+'' 


iV r F ^li 'J /'I — 'J 
^^J,r ^L,-:i/z 


1033.119 


290.2 


49.6 


1.81 


-2) 


2 




^^J,t ^l.d/Z 


1892.227 


158.4 


140.4 


5.9i 


-2) 


2 




/V T^' Tj' T n /I 1 -I -I •) / '1 

^*KaKf,,J ^U2,o/2 -^ll.d/2 


746.3 


401.7 


89.3 


5.5( 


-4) 


3 




^Ka,Ki,,J =lll,3/2-OoO,l/2 


1115.204 


268.8 


53.5 


3.1( 


-2) 


6 




NKaK^,J =lll,l/2-000,l/2 


1139.6 


263.1 


54.7 


2.9{ 


-2) 


3 




NKaK^,J =3l2,5/2-3o3,5/2 


999.8 


299.8 


223.9 


2.3( 


-2) 


3 


H3O+ 


■/jf,+=Oo,--lo,+ 


984.712 


304.4 


54.7 


2.3( 


-2) 


1 




>^if, + =43,+-33,- 


1031.294 


290.7 


232.2 


5.1( 


-3) 


1 




•^X, + =42,+-32,- 


1069.827 


280.2 


268.8 


9.8( 


-3) 


1 




Jk, + =&2,--&2,+ 


1454.563 


206.1 


692.6 


7.1( 


-3) 


1 




•^i<:,+=2i,--2i,+ 


1632.091 


183.7 


143.1 


1.7( 


-2) 


1 


CO 


10-9 


1151.985 


260.2 


304.2 


1.0( 


-4) 


1 




16-15 


1841.346 


162.8 


751.7 


4.1( 


-4) 


1 


"CO 


5-4 


550.926 


544.2 


79.3 


1.1( 


-5) 


1 




10-9 


1101.350 


272.2 


290.8 


8.8( 


-5) 


1 


CI8O 


5-4 


548.831 


546.2 


79.0 


1.1( 


-5) 


1 




9-8 


987.560 


303.6 


237.0 


6.4( 


-5) 


1 




10-9 


1097.163 


273.2 


289.7 


8.8( 


-5) 


1 


HCO+ 


6-5 


535.062 


560.3 


89.9 


1.3( 


-2) 


1 


CH^ 


•^F,P=3/22,--l/2i__|_ 


536.761 


558.5 


25.8 


6.4( 


-4) 


2 




Jf,p=5/23,+-3/22,- 


1661.107 


180.5 


105.5 


3.8( 


-2) 


2 




Jj.,P=5/23,_-3/22,+ 


1656.956 


180.9 


105.2 


3.7( 


-2) 


2 


CH+ 


1-0 


835.138 


359.0 


40.1 


6.4( 


-2) 


2 




2-1 


1669.281 


179.6 


120.2 


6.1( 


-2) 


2 


HCN 


11-10 


974.487 


316.4 


280.7 


4.6( 


-2) 


1 




Nj,F^,F =l2,5/2,7/2-0l,3/2,5/2 


974.478 


307.6 


46.8 


6.9( 


-3) 


2 
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Table 2 — Continued 



Species 


Transition 


1/ 


A 




A 


Ref. 






(GHz) 


(/im) 


(K) 


(s 








Nj,Fi,F =ll,3/2,5/2-0l,3/2,5/2 


999.973 


299.8 


48.0 


5.2 


-3) 


2 


NH+^ 


Jp =3/2_-l/2+ 


1012.540 


296.1 


48.6 


5.4 


-2) 


7 




Jp =3/2+-l/2_ 


1019.211 


294.1 


48.9 


5.5 


-2) 


7 


NH3 


lo-Oo 


572.498 


523.7 


27.5 


1.6 


-3) 


1 


C+ 


^P3/2-^Pl/2 


1900.537 


157.7 


91.3 


2.3 


-6) 


5 




Nj,F = 23_7/2-l2,5/2 


1082.909 


276.8 


77.2 


9.8 


-3) 


2 




^J,F =34,9/2-23^5/2 


1632.518 


183.6 


155.6 


3.6 


-2) 


2 




^J,F =l2,6/2-0l,3/2 


526.048 


569.9 


25.2 


9.6 


-4) 


2 


CS 


11-10 


538.689 


556.5 


155.1 


3.3 


-3) 


1 


SH^ 


Ar;^ =3+1-2-1 


1447.012 


207.2 


640.6 


8.1 


-3) 


1 


H 


H{20)a 


764.230 


401.7 











PACS lines 



H20 




221- 


2l2 


1661.008 


180.5 


194.1 


3.1 


-2) 


1 






2l2- 


loi 


1669.905 


179.5 


114.4 


5.6 


-2) 


1 






3o3- 


2l2 


1716.770 


174.6 


196.8 


5.1 


-2) 


1 






3l3- 


2o2 


2164.132 


138.5 


204.7 


1.3 


-1) 


1 






330- 


321 


2196.346 


136.5 


410.7 


6.6 


-2) 


1 






4o4- 


3l3 


2391.573 


125.4 


319.5 


1.7 


-1) 


1 






4l4- 


3()3 


2640.474 


113.5 


323.5 


2.5 


-1) 


1 






221- 


lio 


2773.977 


108.1 


194.1 


2.6 


-1) 


1 






322- 


2ii 


3331.458 


90.0 


296.8 


3.5 


-1) 


1 






7o7-6i6 


4166.852 


71.9 


843.5 


1.2 


0) 


1 






818- 




4734.296 


63.3 


1070.7 


1.8 


0) 


1 






9o9- 


818 


5276.520 


56.8 


1323.9 


2.5 


0) 


1 


OH 


n,j = 


_ 1 3 
2 ' 2 


1 1 

2 ' 2 


1834.747 


163.4 


269.8 


6.4 


-2) 


1 










1837.817 


163.1 


270.1 


6.4 


-2) 


1 




n,j = 


_ 3 5 
2 ' 2 


3 3 
2 ' 2 


2509.949 


119.4 


120.5 


1.4 


-1) 


1 










2514.317 


119.2 


120.7 


1.4 


-1) 


1 




a, J = 


3 7 
2 ' 2 


3 5 
2 ' 2 


3543.779 


84.6 


290.5 


5.1 


-1) 


1 










3551.185 


84.4 


291.2 


5.2 


-1) 


1 




Cl,J = 


_ 1 1 
2'2 


3 3 
2'2 


3786.170 


79.2 


181.7 


3.5 


-2) 


1 










3789.180 


79.1 


181.9 


3.5 


-2) 


1 


CO 




14 


-13 


1611.794 


186.0 


580.5 


2.7 


-4) 


1 






16 


-15 


1841.346 


162.8 


751.7 


4.1 


-4) 


1 






18 


-17 


2070.616 


144.8 


945.0 


5.7 


-4) 


1 






22 


-21 


2528.172 


118.6 


1397.4 


1.0 


-3) 


1 






23 


-22 


2642.330 


113.5 


1524.2 


1.1 


-3) 


1 






24 


-23 


2756.388 


108.8 


1656.5 


1.3 


-3) 


1 






29 


-28 


3325.005 


90.2 


2.399.8 


2.1 


-3) 


1 






30 


-29 


3438.365 


87.2 


2564.8 


2.3 


-3) 


1 






31 


-30 


3551.592 


84.4 


2735.3 


2.5 


-3) 


1 






32 


-31 


3664.684 


81.8 


2911.2 


2.7 


-3) 


1 
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Table 2 — Continued 



Species 


Transition 


U 

(GHz) 


A 


(K) 


A 
(s-i) 


Ref. 




33-32 


3777.634 


79.4 


3092.5 


3.0(-3) 


1 




36-35 


4115.606 


72.8 


3668.8 


3.6(-3) 


1 


O 


3Pl-3P2 


4166.852 


63.2 


227.7 


8.9(-5) 


5 




3Po-3P2 


5276.519 


145.5 


326.6 


1.8(-5) 


5 


c+ 


■^P3/2^'^Pl/2 


1900.537 


157.7 


91.3 


2.3(-6) 


5 



^Fine/hyperfine-structure transition. Only tlie strongest component is 
listed here. 



Reference s. — (1) JPL spe ctral line catal o gue, Pickett et al.l ll l998h: 
(2) CD MS, iMiiller et al] l l200lf) : (3) iBrudereJ i200e ): (4') ISruderer et al.l 



I 2010bl') : (5) NIST atomic spectra database; (6) iMurtz et al.l lll998l') : (1) 



iHiibers et all l l200g^ 
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Given a temperature and density structure, the molecular excitation and radiative trans- 
fer in the line is computed at each position in the envelope. The resulting sky brightness 
distribution is convolved with the beam profile. A trial abundance of water is chosen (see ex- 
ample in Fig. and is adjusted until the best agreement with observational data is reached. 
The velocity structure is represented either by a turbulent broadening width that is constant 
with position or by some function, for example an infall velocity profile. 

The alternative, 'forward modeling' approach starts from a full physico-chemical model 
and computes the water emission at different times in the evolution for comparison with 
observations (see Fig. 1 of Doty et al. 2004 for a flow-chart of both procedures). In this case, 
one obtains best fit model parameters such as the timescale since evaporation (ofte n labeled 
as the 'age' of the source) or the cosmic-ray ionization rate (e.g.. lDoty et al.ll2006[ ). 



The spherically symmetric models are an important first step but initial Herschel results 
show that they are generally not sufficient to interpret the data. Other c omponents such as 



shocks or UV-heated cav ity walls need to be added in a 2D geometry (Ivan Kempen et al. 
2010bl . IVisser et al.ll201l[ l. 



3. Observations 



3.1. Source and line selection 



The WISH program contains about 80 sources covering a range of luminosities and 
evolutionary stages and uses about 425 hr of Herschel time. The sources are summarized in 
Table [H where they are subdivided in a set of subprograms. Multi-line pointed observations 
at the source position are performed using the HIFI and PACS instruments. In addition, 
small maps over a few arcmin region are made for selected sources using either the on-the- 
fiy (HIFI) or raster mapping (PACS) strategies. The number of sources per (sub)category 
ranges from two for the cold line-poor sources to more than ten for warm line-rich objects, 
large enough to allow individual source pecularities to be distinguished from general trends. 
These deep integrations and thorough coverage of the various types of H2O lines will set 
the stage to design future Herschel programs of larger, more statistically significant samples 
using fewer lines. Specific source selection for each subprogram is discussed in §4. More than 
90% of our sources are visible with ALMA {6 <40°). 

The selection of H2O lines observed with HIFI and PACS is summarized in Fig. [T]and 
Table [2] and is based on information from exis ting ISO-LWS, SWAS and Odin data, and on ex- 
tensive modeling performe d for shocks (e.g., 



19961 . iGiannini et al.ll2006l ). low- mass fICeccarelli et al 



Neufeld fc Dalgarno 



2000a 



1989 



Kaufman fc Neufeld 



van Kempen et al.ll2008l ) and 
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high-mass YSO envelopes (jPoty fc Neufeldlll997l . lWalmsley fc van der Takll2005l . |Poelman fc van der Tak 



20071 ). Three sets of hnes can be distinguished. The ground-state o— and p— H2O hues at 
557 and 1113 GHz are the prime diagnostics of the cold gas. These and other lines connected 
with ground-state levels (e.g., 2i2 — loi at 1670 GHz) usually show strong self-absorption or 
they can even be purely in absorption against the strong continuum provided by the source 
itself. The second group are "medium- J" lines originating from levels around 100-250 K 
above ground, which probe the warm gas (e.g., 2o2 — In at 988 GHz and 812 — 893 at 1097 
GHz). The third group are the highly excited lines originating from levels > 300 K, which 
are only populated in strong shocks (e.g., higher- lying backbone lines) or which are anoma- 
lously excited by collisional or infrared pumping leading to maser emission. Because of high 
optical depths, optically thin isotopic lines are crucial for the interpretation so that deep 
integrations on a number of Hg^O and H^^O lines are included in the program. 

For all deeply embedded YSOs, a common strategy is adopted by observing the ground- 
state o— H2O liQ — loi 557 GHz andp— H2O In — Oqo 1113 GHz lines, as well as the medium- J 
J9-H2O 2o2 - 111 at 988 GHz, 2u - 2o2 at 752 GHz, and the 0-H2O 812 - 221 1097 GHz 
lines, together with at least one Hg^O, OH and CO line (see Table [2]). Several lines of ^^CO, 
C^^O and hydrides are covered serendipitously in the same settings. Higher-J H2O, OH, 
[O I] and CO lines are probed with a series of PACS settings on all low- and intermediate- 
mass sources, which provide 5x5 mini-maps on a 47" x 47" scale with 9.4" spacing to 
probe extended emission. The PACS data also give information on the dust continuum 
emission at wavelengths poorly sampled to date. For high-mass YSOs and a few bright 
low- mass YSOs, full PACS spectral scans are taken, which cover many high-excitation H2O 
lines in an unbiased way. Integration times are generally such that H2O abundances down 
to ~ 10~^ — 10~^° are probed, so that regions with significant freeze-out can be studied. 
Typically 5-6 hr are spent per source to cover all lines with integration times per line ranging 
from 10 minutes to 1 hr. For the pre-stellar cores and protoplanetary disks, long integration 
times with HIFI of up to 12 hr per line are chosen to ensure that any non-detections are 
significant. These are the deepest observations carried out with HIFI across all KPs. 



3.2. Observations and data reduction 

3.2.1. HIFI 

The bulk of the HIFI data are taken in dual-beam switch mode with a nod of 3' using 
fast chopping. The HIFI receivers are double sideband with a sideband ratio close to unity. 
For line-rich sources (mostly high-mass YSOs), the local oscillator was shifted slightly for 
half of the integration time to disentangle lines from the upper and lower side bands. Two 
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polarizations, H and V, are measured simultaneously and are generally averaged together to 
improve the signal-to-noise. In some cases, differences of order 30% are found between the 
two polarizations in which case only the higher quality H-band spectra are used for analysis 
since the mixers have been optimized for H band. Two back-ends are employed: the low- 
resolution back-end (WBS) with an instantaneous bandwidth of 4 GHz at 1.1 MHz spectral 
resolution (~1100 and 0.3 km s~^ at 1 THz, respectively) and the high-resolution back-end 
(HRS) with a variable bandwidth and resolution (typically 230 MHz bandwidth and 250 kHz 
resolution, or ~63 and 0.07 km s~^ at 1 THz, respectively). 

Data reduction of the HIFI spectra involves the usual steps of checking individual expo- 
sures for bad spectra, summing exposures, taking out any basehne ripples, fitting a low-order 
polynomial baseline, and making Gaussian fits or integrating line intensities as appropriate. 
The data are reduced within the Herschel Interactive Processing Environment (HIPE) (Ott 
2010) and can be exported to CLAS^ after level 2 for further analysis. The main beam 
efficiency has been determined to be around 0.76 virtually independent of frequency, except 



for a 15% lower value around 1.1-1.2 THz (jOlbergI l2010l ). and the absolute calibration is 



currently estimated to be better than ~15% for HIFI Bands 1, 2 and 5, and ~30% for Bands 
3, 4, 6 and 7. The rms noise in the WBS data is generally lower than that in the HRS data 
by a factor of 1.4 when binned to the same spectral resolution, due to a a/2 loss factor in 
the HRS autocorrelator. 



3.2.2. PACS 

PACS is a 5x5 array of 9.4" x 9.4" spaxels (spatial pixels) with very small gaps between 
the pixels. Each spaxel covers the 53-210 /xm wavelength range with a spectral resolving 
power ranging from 1000 to 4000 (the latter only below 63 /xm) in spectroscopy mode. In 
one exposure, a wavelength segment is observed in the first order (105-210 /im) and at 
the same time in the second (72-105 yum) or third order (53-72 yum). Two different nod 
positions, located 6' in opposite directions from the target, are used to correct for telescopic 
background. Data are reduced within HIPE. The uncertainty in absolute and relative fluxes 
is estimated to be about 10-20%, based on comparison with the ISO-LWS data. 

The diffraction- limited beam is smaller than a spaxel of 9.4" at wavelengths shortward 
of 110 /im. At longer wavelengths, the point-spread function (PSF) becomes significantly 
larger such that at 200 /im only 40% of the light of a well centered point source falls on 
the central pixel. Even at 100 yum, 30% of the light still falls outside the central spaxel. 



^http : //www ■ iram . f r/I RAM/GI LDAS| 
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Thus, observed fluxes reported for a single pixel have to be corrected for the point-source 
PSF using values provided by the Herschel Science Center. For extended sources, fluxes can 
be summed over the spaxels to obtain the total flux. In case of a bright central source with 
extended emission, such as along an outflow, fluxes at the outflow positions were corrected 
for the leaking of light from the central spaxel into adjacent spaxels. 



3.3. Archival value 



All reduced HIFI and PACS spectra will be delivered to the Herschel Science Center 
approximately one year after the last KP data have been taken, assuming that reliable 
calibration is available in a timely manner. The delivery will consist of reduced and calibrated 
HIFI spectra, line profile parameters (either Gaussian fits or integrated intensities, depending 
on line) and reduced small maps, most likely in FITS format. A quick-look spectral browser 
will be on the team website. For PACS, a data cube of calibrated images will be delivered, 
similar to that for integral field spectrometers. 

Over the past decade, the WISH team has collected a wealth of complementary data on 
various molecular lines for the majority of our sources using the JCMT, IRAM 30m, APEX, 
CSO, Onsala and other single-dish telescopes. These data include tracers of the cloud velocity 
structure and column density {^^CO, ^^CO, C^^O J=l-0, 2-1, and 3-2, probing ~ 30 K 
gas), high temperature lines (CO and isotopologues 6-5, 7-6, probing ~50-100 K gas) and 
high density tracers (CS and C^^S 2-1, 5-4, 7-6; HCO+ and Hi3C0+ 3-2, 4-3, probing 
n > 10^ cm~^), as well as the atomic [C I] line (probing UV radiation). Moreover, grain 
surface chemistry products are observed through various H2CO and CH3OH lines, whereas 
shocks powerful er iough to disrupt g r ain cor e s are traced through S i O emission. Example s 



of such studies are van der Tak et al. 



Nisini et all (120071 ). Marseille et al.l ( 120081 ) 



(l2000bh . lJe)rgensen et al.l (l2004f ). lFuente et all ( l2005b[ ) 



A recent development is the use of heterodyne array receivers such as the 9 pixel HERA 
230 GHz array at the IRAM 30m; the HARP-B 16 pi xel 345 GHz receiver at the JCMT 



[e.g. 



van Kempen et al.l l2009bl . iKristensen et al.l l2010al): and the dual-frequ ency 2x7 pixel 
CHAMP+ 650/850 GHz receiver at APEX (e.g.. Ivan Kempen et al.lboOQdI fl). The spatial 



resolution of these telescopes is typically 7-15", i.e., comparable to that of Herschel at the 
higher frequencies. These data also provide an 'image' of the warm dense gas contained 
within the ~ 37" Herschel beam at 557 GHz. Millimeter interferometry data down to ~1" 
resol ution also exist for a significant subset of the sources in various molecules and lines 



Fuente et al.l l2005aL Ivan der Tak et al.l l2006l . iBenz et al.ll2007l . Ijcargensen et al.l 12009 



(e.g., 

and references cited). Except for the compact hot core, however, such data cannot be 
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compared directly with single dish data from Herschel since a significant fraction of the 
extended emission is resolved out in the interferometer. 

Finally, submillimeter continuum maps obtained with JCMT-SCUBA, IRAM-MAMBO, 



(e.g., 


Di Francesco et al. 


2008 


), as are e 


(Gibb et al. 


2nn4b. 


Boogert et al. 


2008) 



for the three outfiow targets studied in detail (INeufeld et al.ll2009l ). 



The ground-based complementary data will be summarized in a table at the KP Web 
site. For several subprograms, the most relevant data will be available directly at the KP 
web site. Specifically, in collaboration with the 'Dust, Ice and Gas in Time' (DIGIT) KP, 
a searchable web-based archive of complementary submillimeter single-dish, interferometer 
and infrared spectra is being built for low-mass YSOs. 



4. Results 

In the following, a brief description of the scientific motivation and observing strategy 
for each of the subprograms is given together with initial results. Full source and line lists 
are given in Tables [1] and O Unless otherwise specified the cited noise limits are for 0.2-0.4 
km s~^ velocity bins and all stated results refer to main beam antenna temperatures, Tmb- 



4.1. Pre-stellar cores 



Motivation. In recent years, a number of cold, highly extincted cloud s have been iden- 
tified which have a clear central density condensation (see overview by iBergin fc Tafalla 
20071 ). Several of these so-called pre-stellar cores are thought to be on t he verge of collapse 



and thus represent the earliest stage i n the star- formation process (e.g., iTafalla et al.lll998 



Caselli et al.l |2002| . ICrapsi et al.l l2005l ) . It is now widely ac cepted that most molecules are 



highly depleted in the inn er dense parts of these cores (e.g., iCaselli et al.lll999l . iBergin et al. 
2OO2I . iTafalla et al.l l2006l ) . H2O should be no exception, although its formation as an ice 
already starts in the general molecular cloud phase. The physical properties of these cores 
are very well constrained from ground-based continuum and line data and they are the most 
pristine types of clouds, without any internal heat sources and almost no turbulence. Thus, 
they form the best laboratories to test chemical models. Since H2O is the dominant ice 
component and has a very different binding energy and formation route than CO, HIFI 
observations of H2O will form a critical benchmark to test gas-grain processes. 
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Fig. 4. — HIFI spect ra of the o— H2O Im — loi 557 GHz line. From top to bottom: the 
pre-stellar core B68 (ICaseUi et al.ll2010l ). the low-mass Sta ge YSO L1527, th e low- mass 
Stage I YSO TMR-1, and the protoplanetary disk DM Tau JSergin et al.lboiObh . 
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HIFI will be key to address the questions: (i) Where does the onset of H2O ice for- 
mation and freeze-out occur in dense pre-stellar cores? (ii) How effective are non-thermal 
desorption mechanisms (photodesorption at edge, cosmic-ray induced desorption in center) 
in maintaining some fraction of O and H2O in the gas phase? (iii) Is the H2O abundance in 
the gas phase sensitive to the pre-stellar core environment? 

Source and line selection: Two well-studied pre-stellar cores, B68 and L1544, are ob- 
served in the o— H2O 557 GHz line at the central position down to ~2 mK rms. Both are at 
d <200 pc and have bright millimeter continuum dust emission, a simple morphology and 
centrally concentrated density profiles. B68 is an isolated dark globule exposed to the general 
interstellar radiation field whereas L1544 is a dense core embedded in a larger cloud. The 
sensitivity lirnits ar e based on H2O column densities estimated using the model results of 



Aikawa et al.l ( 120051 ) for chemically 'old' sources like L1544. From the predicted abundance 
profile and the density structure, column densities have been determined at different impact 
parameters and convolved with the HIFI beam at 557 GHz. T he brightness temperatures 



(2000 


) and 


Keto et al. 


(2004) 



of the core, the initial H2O abundance in the chemical models, as well as the assumed 
ortho/para ratio of both H2O and H2. Pie-Herschel model p redictions ranged from peak 



antenna temperatures of below 10 mK ( iHoUenbach et al.ll2009[ ) to several hundred mK 



Initial results. ICaselli et al.l (|2010[ ) present early results for the two cores, with the B68 
0-H2O lio — loi spectrum included in Fig. HJ These spectra took ~4 hr integration (on-l-off ) 
each. No emission line is detected in B68 down to 2.0 mK rms in 0.6 km s~^ bins. This is 
more than an order of magnitude lo wer than the previous upper limit obtained by SWAS 
for this cloud (IBergin fc Snelll 120021 ) as well as most pre-Herschel model predictions. The 
3a upper limit corresponds to an o— H2O column density < 2.5 x 10^^ cm~^ and a mean 
line-of-sight abundance < 1.3 x 10~^. This is consistent with the non detection predicted by 
the model of iHoUenbach et al.l ( 120091 ) applied to B68. 



Interestingly, the L1544 spectrum shows a tentative absorption feature against the sub- 
millimeter continuum, which closely matches the velocity range spanned by the CO 1-0 
emission line. If the water absorption is confirmed by deeper integrations which are currently 
planned within WISH, it provides a very powerful tool to determine the water abundance 
profile along the line of sight. 



Initial radiative transfer analyses, based on the iKeto fc Casellil ( I2OIOI ) models, imply an 
H2O abundance profile for both clouds which peaks at an abundance of ~ 10~^ in the outer 
region of the core, where the density is too low for significant freeze-out but where the gas is 
sufficiently shielded from photodissociation [Ay > 2 mag). The abundance drops by more 
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Fig. 5. — Cartoon of a protostellar envelope with the different physical components and 
their nomenclature indicated. IWS s tands for internal working surfaces. The indicated scale 
is appropriate for a low- mass YSO (IVisser et al.ll201ll ). On this scale, the ~100 AU radius 
disk surrounding the protostar is not visible. 
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than an order of magnitude in the central part due to freeze-out. These new predictions are 
based on an improved physical structure of the L1544 core and on an H^O abundance p rofile 
which includes the freeze-out, a generalized desorption from [Roberts &: Millarl (120071 ) and 
photodissociation, but no detailed gas phase and su rface chemistry. Suc h abundance profiles 
differ from those expected from chemical models by lAikawa et al.l ( l2005l ). who predict higher 
H2O abundances toward the core center and thus brighter lines. Investigations are currently 
under way to understand which processes are producing the different abundance structures. 

These results demonstrate that questions (i) and (ii) can be addressed by the Herschel 
observations, provided the integrations are deep enough. Question (iii) requires observations 
of a larger sample of cores besides those studied in WISH. 



4.2. Low-mass young stellar objects 



Motivation. Once the cloud has started to collapse and a protostellar object has formed, 
the envelope is heated from the inside by the accretion luminosity. The physical structure 
of these low-mass deeply embedded objects 100 Lq) is complex, with large gradients 
in temperature and density in the envelope, a compact disk at the center, and stellar jets 
and winds creating a cavity in the envelope and affecting the surroun ding material (e.g.. 



van Dishoeck et al.l Il995l |j0rgensen et al. 



2005al IChandler et al.l 120051 ) (see Fig. E]). The 



high spatial and spectral resolution at far-infrared wavelengths offered by Herschel are key 
to disentangling these processes. 

Several evolutionary stages can be distinguished for low-mass YSOs (see nomenclature 
by lRobitaille et al.ll2006l ). In the earhest deeply embedded Stage 0, the main accretion phase, 
the envelope mass is still much larger than that of the protostar or disk {M^m » ^* and 
>> Afdisk)- At the subsequent Stage 1, the star is largely assembled and now has a much 
larger mass than the disk (M^, >> Mjisk), but the envelope still dominates the circumstellar 
material (Mgnv > Miisk)- At Stage 2, the envelope has completely dissipated and only a 
gas-rich disk is left around the young pre-main sequence star. Observationally, low-mass 
YSOs are usually classified based on their SEDs according to their spectral slo pes (as in the 
Lada classes) or bolometric temperatures (see summary in lEvans et al.ll2009l ). The Class 
and II sources generally correspond well to the evolutionary Stage and 2 cases. In contrast, 
the traditional Class I sources are found to be a mix of true Stage 1 sources and of Stage 2 
disks which are either seen edge-on or which have a lar ge amount of foreg round extinction 
resulting in a rising infrared spectrum (see discussion in lCrapsi et al.ll2008l ). 



ISO-LWS observations have shown that low-mass YSOs are copious water emitters, with 
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Fig. 6.— HIFI spectra of the p-HaO In - Oqo 1113 GHz (left) and 2o2 - In 988 GHz 
fright ) lines. From top to b ottom: the high-mass YSO W3 IRS5 (L = 1.5 x 10^ Lq, d =2.0 



kpc) (IChavarrfa et a 



pc) (IJohnstone et al 



2ninh . the intermediate-mass YSO NGC 7129 FIRS2 (430 Lq, 1260 
gOloh . and the low-mass YSO NGC 1333 IRAS2A (20 Lq, 235 pc) 



(IKristensen et al.ll2010bl ). All spectra have been shifted to a central velocity of km s^^. 
The red-shifted absorption features seen in the In — Oqo spectrum toward W3 IRS5 are due to 
water in foreground clouds whereas the small emission line on the blue side can be ascribed 
to the SO2 139 5 — 128 4 transition. 
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more than 70% of the St age sources dete cted in H2O hnes with energies up to 500 K above 
ground (see summary by lNisini et aLll2002l ). Two scenarios have been put forward to explain 
these strong hnes. In one model, the lines arise primarily from the bulk of the envelope, with 
the lower-energy li nes originating in the outer envelope and the higher-energy lines from 
the inner hot core (IMaret et al.ll2002l ). In the alternative interpretation, most of the H2O 
emission is associate d with the outflows impacting the surrounding gas on larger scales (e.g., 
Giannini et al.ll200ll ). 



In contrast, none of t he low-mass Stage 1 sources show H2O lines in the ISO-LWS 



beam (iGiannini et al.ll200l[ ). Possible explanations given in that paper include dissociation 
resulting from enhanced penetration of X-rays or UV photons or a longer timescale enhancing 
freeze-out. HIFI can distinguish between the different scenarios by providing much deeper 
searches than ISO-LWS, as well as the photodissociation products OH and O, and X-ray 
products such as ions. Indeed, OH/H2O abundance ratios are expected to change by factors 
of >100 due to UV and X-rays (IStauber et al.ll2006l ). Complementary infrared data on H2O 
ice are available for many of these sources so that the gas/ice ratio can be measured directly. 

Key questions to be addressed are: (i) What is the origin of the strong H2O emission 
in low-mass Stage protostars? Passively heated envelope or shocked outflow material? 
What is the role of the outflow cavity? (ii) Is the high water abundance governed by high- 
temperature chemistry or grain mantle evaporation? (iii) Do low-mass Stage 1 sources 
indeed have much lower H2O abundances? If so, is this due to enhanced enhanced UV, 
X-rays or other effects such as lower outflow activity? (iv) What is the HDO/H2O ratio 
in the protostellar phase, and what does this tell us about the efficiency of grain surface 
deuteration? How does HDO/H2O evolve during star formation? Is it similar to that found 
in comets? 

Source and line s election. A set of 14 deeply embedded Stage sources has been se- 
lected from the list o f lAndre et al.l (I20001) and 15 Stage 1 sources selected from the lists of 
Tamura et al.l (Il99ll ). lAndre fc Montmerld (119941) and the spectroscopic part of the Spitzer 
'Cores to Disks' legacy program ( lEvans et al.l l2003l : |2009[ ). All sources have d < 450 pc 
and are well ch aracterized with a var i ety o f other single-dish submillimeter telescopes and 
interferometers (|j0rgensen et al.ll2004j : 120071 ) . In particular, only bonafide Stage 1 embedded 
YSOs are inc luded, with the stage verified through compleme ntary single-dish and interfer- 
ometer data (Ivan Kempen et al.ll2009bl . |j0rgensen et al.ll2009l ). 



The line selection for the Stage sources follows the general philosophy outlined in § 13. H 
with typical integration times for the HIFI lines such that 20-100 mK rms is reached in ~0.5 
km s~^ bins depending on the line, and down to 5 mK for Hg'^O. For reference, observed 
widths of molecular lines originating in the bulk of the envelopes are typically 1-2 km s~^ 
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FWHM. The noise limits are based on ID envelope models by Ivan Kempen et al.l (120081 ) 
for L483, one of the weaker Stage sources, and do not include any outflow contribution. 
Not all lines are observed for all sources: some time-consuming higher frequency lines are 
dropped for the weaker ones. Selected PACS lines are observed to cover the higher excitation 
emission and obtain information on the spatial extent. Full PACS spectral scans are made 
of four Stage sources (NGC 1333 IRAS 4A/B, IRAS2A, and Serpens SMMl) to obtain an 
unbiased view of the higher H2O lines. For the Stage 1 objects, only four H2O lines and one 
CO line are observed with HIFI. The set of PACS lines observed for these sources is also 
more limited. 

Initial results. Figure [6] presents spectra of the lowest two P-H2O lines toward one of the 
brightest low- mass Stage sources, NGC 1333 IRAS2A flKristensen et al.ll2010bl ). whereas 
Fig. m includes the lowest 0-H2O line toward the Stage source L1527. The profiles are 
surprisingly broad and complex. They can be decomposed into three components: a broad 
(>20 km s~^), medium-broad (~5-10 km s~^) and narrow (2-3 km s~^) component, with 
the latter seen primarily in (self-)absorption in the ground-state transitions. The central 
velocities of the three components can vary relative to each other from source-to-source, 
probably due to different viewing angles and geometries of the systems. For one Stage 
source, NGC 1333 IRAS4A, the 2o2 — In line shows an inverse P-Cygni profile indicative of 
large scale infall. Interestingly, the H2^0 lines only reveal the broad component. 

From the profiles themselves, it can be concluded immediately that the bulk of the water 
emission arises from shocked gas and not in the quiescent envelope: a hot core region with 
an enhanced H2O abundance an d an infall velocity profile gives lines with widths of at most 
4-5 km s~^ (jVisser et al.l 120111 ). The broad component is associated with shocks on large 
scales along the outflow cavity walls, as traditionally probed with the standard high velocity 
CO outflows (labelled as 'shell shocks' in Fig. E]). The medium component is linked to shocks 
on smaller scales in the dense inner envelope (< 1000 AU), based on comp arison of the line 
profile s with those of high density tracers and with interferometer data ( iKristensen et al. 
2010b|). The H2O/CO abundance ratios in the shocks are high, ~0.1-1, corresponding to 
H2O abundances of ~ 10~^ — 10~^ with respect to H2. The highest abundances are found at 
the highest velocities, indicating that only a few % of the gas is processed at high enough 
temperatu res where all O is dri ven into H2O. This small percentage is consistent with the 
findings of iFranklin et al.l (120081 ) based on SWAS data. 



The narrow water absorption lines probe the quiescent outer envelope and any fore- 
ground material. Together with the absence of narrow H2^0 lio — loi emission and the 
maximum narrow emission that can be hidden in the H2O profiles, these lines constrain the 
H2O abundance in the bulk of the cold envelope to ~ 10^^. The lack of (narrow) 112*^0 
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Fig. 7.— PACS 5x5 spectral maps of the [O I] ^Pi Pa 63 /xm (left), H2O 2i2 - loi 179 
fim (middle) and CO J=16-15 163 /xm lines (right) toward the low-mass Stage 1 YSO IRAS 
15398-3359. All lines are e xtended along the outflow direction seen in the low- J CO map of 
van Kempen et al.l (j2009al ). 
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2o2 — 111 emission down to ^ 10 mK rms l imits the water abundance in the inner ice evap- 
oration zone to < 10~^ (jVisser et aLll201l[ ). In contrast, the warm quiescent envelope with 
T > 25 K is clearly revealed by HIFI in narr ow C^^O 9-8 emi ssion lines, where the CO 
abundance jumps by a factor of at least a few (lYildiz et al.ll2010l ). 



Figure m includes the 0-H2O 557 GHz lio — loi HIFI spectrum of the Stage 1 source 
TMR-1. This is the first detection of the H2O 557 GHz line in this Stage. The fact that the 
line is weaker than for the Stage source is consistent with the non-detection by ISO. Again, 
the broad line points to an origin in the outflow shocks rather than the quiescent envelope. 

The combination of ground-based, HIFI and PACS data allows the complete CO lad- 
der from J=2-l up to 43-42 to be characterized for a low-mass YSO for the first time 



(Ivan Kempen et al. 



2010bl ). Because CO is much less sensitive to chemical effects than H2O, 



this provides an opportunity to test various gas heating processes. Figure M shows the CO 
excitation diagram for the Stage 1 source HH 46, giving a rotational temperature of 410 ± 25 
K. Also included in this figure are data for another WISH source. Serpens SMMl, for which 
a full PACS spectral scan has been obtained. The entire CO ladder can be fit by a single 
rotational temperature of 520 ± 20 K. CO ladders have also been published for two sources 
from the DIGIT first results, the Stage 1-2 transitional source DK Cha and the disk around 
the Herbig star HD 100546 (Ivan Kempen et al.l l2010al . ISturm et al.ll2010l ). In both cases, 
the data indicate excitation temperatures of several hundred K, but with evidence for a two 
temperature component structure in at least one source. Such a two component fit does not 
necessarily imply two different physical regimes or a temperature gradient, however. Optical 
depth effects in the CO lines can also lead to curvature in the rotational diagrams. 

Because of the limited information that can be derived from an excitation diagram, a 
detailed model has been develope d for HH 46 using the know n physical structure of the quies- 
cent envelope as a starting point (Ivan Kempen et al.ll2009dl ). The passively heated envelope 
can explain only the lowest CO rotational transitions, so additional heating mechanisms are 
needed (see Fig. [5]). One possibility is UV photon heating, in which the UV photons pro- 
duced in the immediate protostellar environment can escape through the cavity carved by 
the jet and wind from the young star and impinge on the envelope at much larger distances 
from the star (up to a few thousand AU). It is found that this UV heating of the outflow 
cavity walls can reproduce well the intermediate transitions up to J < 20, depending on the 
precise PDR model used. Another mechanism is necessary to populate the highest CO levels 
J > 20: small-scale shocks with velocities ~ 20 km s~^ paste d along the walls can explain 



these lines well (Ivan Kempen et al.ll2010bl . IVisser et al.l 120111 ). This model solution is not 
unique and is sensitive to model details, but it demonstrates that a good fit can be found 
with just two free parameters: the UV flux and the shock velocity. 
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Figure [7] presents maps of the CO, H2O and [O I] lines observed with PACS toward 
the isolated Stage 1 source IRAS 15398-3359, known for its peculiar carbon chemistry 
(ISakai et al.l 120091 ) . All lines are clearly ex tended in the outflow direction as mapped in 
lower- J CO lines (Ivan Kempen et al.ll2009al ). confirming their association with the shocks 
and outflow cavities. 

The amount of cooling produced by each of the species can be quantified from the PACS 
data. In the central HH 46 spaxel, [O I] and OH dominate the cooling (at least 60% of the 
total cooling in far-infrared lines) followed by H2O (~ 30%) and CO (~ 10%). The total 
far-infrar ed line cooling of a t least 2.4 x 10~^ L© is consistent with the kinetic luminosity 
fo und bv Ivan Kerapen et al.l (l2009ctf ) from CO outflow maps and follows the relation found 
by Giannini et al. ( 200l[ ) between these two quantities. The [O I] 63 /im line is spectrally 
resolved with PACS, and reveals a jet component with velocities up to 170 km s~^ associated 
with the blue- and red-shifted outflow seen in optical atomic lines in less extincted regions of 
the flow. This is the first time that such a jet has been seen in a far-infrared line, providing 
a powerful diagnostic of the "hidden jet" in the densest optically opaque regions. OH is de- 
tected strongly on source but not along the outflow axis, in contrast with other species. Both 
the [O I] and OH emission are ascribed to a dissociative shock caused by the jet impacting 
on the inner dense (> 10^ cm~^) envelope; disk emission is exclu ded by scaling the observe d 



disk emission from Stage 2 sources to the distance of HH 46 (Ivan Kempen et al.l l2010bl ). 
In dissociative shocks, OH is one of the first molecules to reform after molecular dissocia- 
tion, whereas the UV radiation from the shock and the high H /H^ ratio prevent effective 
H2O formation in the postshock gas (iNeufeld &: Dalgarnolll989l ). Once H2 is reformed, the 
heat of H2 formation mainta ins a temperature plateau i n the warm postshock gas in which 
significant H2O is produced ( iHoUenbach fc McKeelll989l ). 



The OH chemistry has been further investigated by lWampfler et al.l (120 10[ ). The 163 fim 
line seen by PACS toward HH 46 was observed but not detected with HIFI. The upper limit 
implies a line width of at least 11 km s~^. Inspection of the PACS data for six sources reveals 
a correlation between the OH and [O I] fluxes, which, together with the broad profiles, is 
consistent with production in (dissociative) shocks. The relative strengths of the various OH 
far-infrared lines are remarkably constant from source to source, indicating similar excitation 
conditions. Only the 119 /xm lines connecting with the ground state vary because they are 
affected by (foreground) absorption. Due to the lack of detected quiescent OH (and H2O) 
emission in the HIFI spectra, no abundances in the envelope can yet be determined to test 
the model predictions of enhanced OH/H2O ratios due to X-rays. 

The HDO/H2O ra tio has been determined for the low-mass protostar NGC 1333 IRAS2A 
by iLiu &: et al.l ( I2OIOI ). Five HDO lines arising from levels ranging from 20 K up to 170 K 
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have been observed using ground-based telescopes. The inner (> 100 K) and outer (< 100 
K) HDO abundances with respect to H2 are well constrained to (0.7-1.0) x 10^'' (3cr) and 
(0.1-1.8) X 10~^ (3cr) respectively. If the H2O and isotopologue lines discussed above can 
be associated with the same component as seen in HDO, the inferred limit on the warm 
H2O abundance implies a HDO/H2O abundance higher than ~1% in the inner envelope , 



consistent with the value found in the hot corino of IRAS16293 (~3%, iParise et al.ll2005bl ). 
Although absolute abundances of HDO and H2O jump by more than one order of magnitude 
between the cold and warm envelopes, the HDO/H2O ratio in the outer cold envelope is 
found to have a similar value to that of the inner part, ~ 1 — 18%. 

The detection of HDO lines in other low-mass protostellar sources has proven difficult 
from the ground, even under the best weather conditions. Deep ifersc/ie/ observations of high- 
frequency HDO lines coupled with deep Hg^O data are likely needed to pin down the cold 
HDO/H2O ratios. Alternatively, interferometric observations of the higher excitation HDO 
lines, coupled with similar data of the 3i,s — 22n 203 GHz line (the only H2^0 line readily 



imaged from the ground, |j0rgensen fc van Dishoeckl l2010bl ). ca n constrain the HDO/H2O 



ratio in the warm gas in the innermost regions. Initial results by |j0rgensen fc van Dishoeck 



(j2010al ) for the NGC 1333 IRAS4B system indicate much lower values of < 6 x 10 ^ than 



those cited above for IRAS2A. This large difference is not yet understood. 

Overall, it appears that the answer to the first question is that shocks dominate the 
observed water emission in low-mass YSOs rather than the hot core. Both ice evaporation 
and high temperature chemistry are likely important in producing high water abundances 
in shocks. The answers to the other two questions await more data and analysis. 



4.3. Outflows from low-mass young stellar objects 

Motivation. Deeply embedded protostars drive powerful molecular outflows. These out- 
flows extend over arcminute scales for nearby low-mass objects and can be resolved spatially 
with Herschel. The brightest outflow knots are clearly separated, often by arcminutes, from 
the dense protostellar envelopes targeted in the low-mass subprogram ( §4.21) . Outflows ex- 
hibit a great diversity in their velocity, degree of collimation, opening angle, extent, and 



gas kinetic temperature (IBachiller fc Tafallal Il999l ). raising the question of whether these 



characteristics can be linked to the evolutionary state of the driving sources or whether 
other parameters (e.g., structure of the protostellar envelope, environment) play a role. In 
contrast with ground-based observations of low— J CO transitions, Herschel data probe the 
high-excitation regions of the flow, where the 'hidden jet' interacts with surrounding ma- 
terial, and they directly target one of the main coolants, H2O, thus addressing the shock 
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Fig. 8. — Diagram illustrating the CO excitation. The number of molecules in the upper 
level divided by its statistical weight is plott ed versus rotational quantu m number J„. Left: 
the low-mass Stage 1 YSO HH 46 (based on Ivan Kemp en et al.l (l2010bl )). Right: the more 
luminous embedded source SMMl (Goicoechea et al. in prep.). For the latter source, a 
complete PACS spectral scan has been obtained. For both sources, the fluxes from the 
central spaxel have been used. The dashed lines give the thermal distribution at the indicated 
temperatures. The hatched areas in the HH 46 diagram indicate the different physical 
processes thought to be responsible for the CO emission. 




Fig. 9. — Herschel-PACS image of the H2O 179 fim line toward the Stage protostar L 1157, 
together with (spectrally unresolved) PACS spectra at three positions. The water emission 
is most significant in 'hot spots' along the outflow where the (r otating) jet intera cts directly 
with the cloud as well as toward the protostar itself (based on lNisini et al.ll2010l ). 
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physi cs and energetics as a function of distance from the central driving source ( iNisini et al. 



2000). 



Outflows represent a major perturbation in the chemistry of the core gas. In the 
youngest systems, the SiO, CH3OH and H2O abundances are increased by orders of mag- 
nitude as a result of grain destruction and high-temperature gas-phase chemistry. By de- 
termining the H2O abundances in several spatially and kinematically resolved components 
within a given outflow with HIFI, and by doing this for a sequence of outflows of different 
ages, the H2O cycle -from grain mantle, to gas, to mantle- can be reconstructed. This 
allows a specific puzzle posed by ISO-LWS to be addressed: in most protostellar outflow 
regions, the H2O abundance is an order of magnitude lower than the value of 10~^ expected 
if high-temperature reactions drive all available gas-phase oxygen into H2O or if ice mantles 
evaporate. Ionizing and photodissociative effects of radiation emitted by nearby dissociative 
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Key questions to be addressed in this subprogram are: (i) What is the nature of the flow: 
dense molecular jets or shocked ambient gas along the cavity walls? Do the characteristics 
change with distance to the central driving source and with evolution of the system? What 
is the contribution of water emission to the energy lost in outflows? (ii) What is the physical 
impact of the outflows on their surroundings? The transfer of energy and momentum from 
the flow to the envelope and cloud plays a critical role in controlling future star formation, 
(iii) What is the chemical effect of the outflows? Can shock-generated UV radiation modify 
the physics and chemistry of nearby molecular gas? Why is not all oxygen driven into H2O? 

Source and line selection. The outflow subprogram consists of three parts: a survey 
part, a line excitation part and spatial maps. In the survey part, 26 flows from Stage 
and Stage I low-mass YSOs are observed in two II2O lines — the ground-state o— H2O 557 
GHz line with HIFI down to ~20 mK and the 179 fim p— H2O line with PACS to similar 
limits — in two 'hot spots' (blue and red-shifted outflow lobes), providing a statistical basis 
for the assessment of time evolution effects. The driving sources themselves are covered in 
the low- and intermediate-mass YSO subprograms (see §4.21 and 14. 4p . The line excitation 
survey targets a sub-sample of three Stage outflows for which H2O emission has previously 
been detected (L1157, L1448, VLA1623). These sources are observed in 15 spectral lines 
with HIFI and PACS (H2O, OH, [O I] and CO) to characterize their physical structure and 
energy budget. Integration times are such that typically 20-40 mK rms is reached. In the 
third part, these same three flows are spatially mapped with HIFI down to ~20 mK using 
the on-the-fly (OTF) mode in the H2O 557 GHz line to determine the spatial distribution 
of shocked H2O and compare it with other outflow tracers like CO, SiO and CH3OH. In 
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addition, they are imaged with PACS in the 179 /im hne. For NGC 2071, whi ch is known 
from Spitzer data to be particularly rich in H2O emission (IMelnick et al.ll2008l ). full PACS 
spectral scans are performed at a few positions. 

Initial results. Figure [9] shows the Herschel-PACS image of H2O in the 179 /im line 
towar d the Stage sour ce LI 157 obtained during the science demonstration phase of Her- 
schel (INisini et al.ll2010l ). PACS spectra at selected positions are included, illustrating the 
high S/N of the data. The map beautifully reveals where the shocks deposit energy into the 
molecular cloud, lighting up the water emisson along the outflow. The H2O e mission is spa 



tially correlated with that of pure rotational lines of H2 observed with Spitzer ( iNeufeld et al. 



20091 ). and corresponds well with the peaks of other shock-produced molecules such as SiO 
and NH3 along the outflow. In contrast with these species, H2O is also strong at the source 
position itself. The analysis of the H2O 179 fim emission, combined with existing Odin and 
SWAS data, shows that water originates in warm compact shocked clumps of few arcsec in 
size, where the water abundance is of the order of 10~^, i.e., close to that expected from 
high-temperature chemistry. The total H2O cooling amoun ts to 23% of the tot al far-IR 
energy released in shocks estimated from the ISO-LWS data (IGiannini et al.ll200lh . 



The PACS spectra contained in Figure |9] are spectrally unresolved. HIFI spectra of 
H2O and CO at the Bl position in t he blue outflow lobe (see Fig. [9]) have been presented 
and analyzed by iLefloch et al.l ( I2OIOI ) as part of CHESS, showing th at the H2O abundance 
increases with shock velocity (see also §4.21 and iFranklin et al.l ( 120081 )). ICodella et al.l (l2010l ) 
propose that this reflects the two mechanisms for producing H2O: ice evaporation and high 
temperature chemistry where the O -|- H2 and OH + H2 reactions become significant (see 
§ 12. H and Fig. |3]), with the latter reactions becoming more important at h i gher v elocities. The 
presence of NH3 (another grain surface product) seen by ICodella et al.l (120101 ) only at lower 
velocities is consistent with this picture. The WISH data can test this proposed scenario at 
many more outflow positions. 

In summary, the limited data to date illustrate the potential of water to probe and 
image shock physics and chemistry along the outflows. 



4.4. Intermediate-mass young stellar objects 

Motivation. Intermediate-mass YSOs with ~ 10^ — 10^ Lq (stellar masses of a few 
to 10 M0) form the bridge between the low- and high-mass samples. Besides covering 
a significant luminosity range between the two well-studied extremes, intermediate-mass 
YSOs are interesting in their own right. For example, they are often found in clustered 
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environments, so they can be used as laboratories for investigating high-mass cluster regions 
which are often more confused due to their larger distances. Intermediate mass YSOs may 
be in an interesting regime where heating from t he outside by radiatio n from neighboring 
stars is significant compared with internal heating ( 1j0rgensen et al.ll2006l ). This would result 
in a more extended warm zone and a different chemical structure of the envelope than that 
of their low- and high-mass counterparts. Also, emerging intermediate-mass YSOs have 
associated PDRs or reflection nebulae indicating that not only shocks but also UV photons 
play a role in clearing the surroundings. Finally, they are the prec ursors of Herbig A e/Be 
stars, whose protoplanetary disks have been very well studied (e.g., iMeeus et al.ll200l[ ). but 
for which very little is known about their formation history. 

Key questions are therefore: (i) What is the structure of the envelopes of intermediate- 
mass YSOs? Do they have a different temperature and chemistry structure compared with 
their low- and high-mass counterparts? Does the enhanced UV compared with the low-mass 
counterparts play a role? (ii) How do intermediate mass YSOs evolve? To what extent are 
the answers to the questions in § 14.21 14. 3[ 14.51 different for these sources? 



Source and line selection. Although intermediate mass YSOs are more energetic than 
their low-mass counterparts, they are not bright enough to observe at the distances of their 
high-mass cousins. Surprisingly few nearby (< 400 pc) candidates are known, so our team has 
conducted an active search over the past years to identify a representative sample to observe 
with HIFI. Six intermediate mass YSOs with distances <1 kpc have been selected from 
this sample and a significant amount of complementary submillimeter line an d continuum 
observations have been obtained (e.g.. iFuente et al.ll2005bl . ICrimier et al.ll2010l ). All of these 
sources are still at the deeply embedded stage, comparable to Stage in the low-mass case. 
The Herschel line selection follows the general criteria outlined above and the noise levels 
are similar to those for the low-mass YSOs. 

Initial results. Observ ations of va r ious water lines toward NGC 7129 FIRS2 (430 L©, 
1260 pc) are presented in iFich et al.l (l2010l ). IJohnstone et al.l (120101 ) and the two lowest 
H2O lines detected with HIFI are included in Fig. [61 The line profiles are remarkably 
similar to those of the low-mass YSO NGC 1333 IRAS2A, showing both a broad component 
(FWHM ~25 km s~^) due to shocks along the outflow cavity as well as a medium-broad com- 
ponent (FWHM ~6 km s~^) associated with small-scale shocks in the inner dense envelope. 
Quantitative analysis of all the H2O and isotopologue lines observed for this object result 
in an outer envelope (< 100 K) abundance of water of order ~ 10~^, an order of magnitude 
higher than for the low- and high-mass YSOs. Alternatively, both the high-J CO and H2O 
lines seen with HIFI and PACS can be analyzed in a single slab model with a temperature 
of ~ 1000 K and density ~ 10'' — 10^ cm~'^, representative of a shock in the inner dense 
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envelope (IFich et al.ll2010[ ). This leads to a typical H2O/CO abundance ratio of 0.2-0.5. 

The results for this first source demonstrate that the observing strategy for this sub- 
program is sound and that the data can provide the quantitative information to answer the 
above questions, once more sources have been observed and analyzed. 



4.5. High-mass young stellar objects 

Motivation. High-mass stars (> 10 Mq, > 10^ Lq), though few in numbe r, play a major 



role in the energy budget and shaping of the Galactic environment (e.g., ICesaroni et al. 



20051 ). These are also the type of regions that dominate far-infrared observations of starburst 
galaxies. Their formation is not at all well understood due to large distances, short timescales 
and heavy extinction. 

Recent results show that the embedded phase of the formation of O and B stars can 
be divided empirically into several classes. The earliest stage are the massive pre-stellar 
cores (PSC), which are local temperature minima and density maxima and which represent 
the initial conditions. Their observational characteristics are large column densities, low 
temperatures, and absence of outflow or maser activity. The next stage are the High-Mass 
Protostellar Objects (HMPOs), where the central star is surrounded by a massive envelope 
with a centrally peaked temperature and density distribution. HMPOs show signs of active 
star formation through outflows and/or masers. This subprogram distinguishes mid-IR- 
bright and -quiet objects (mIRb / mIRq), with a boundary of 10 Jy at 12 fj,m. The next 
evolutionary phase is represented by Hot Molecular Cores (HMC), with large masses of warm 
and dense dust and gas, and high abundances of complex organic molecules evaporated off 
dust grains. HMCs have regions of T > 100 K that are >0.1 pc in size. The final stage 
considered here are Ultracompact HII regions (UCHII), where large pockets of ionized gas 
have developed but which are still confined to the star. UCHII sources have free-free flux 
densities of more than a few mJy (for a distance of a few kpc) and sizes of 0.01-0.1 pc. As with 
the previous category, this criterion may imply a luminosity of >10^ Lq. The embedded phase 
ends when the ionized gas expands hydrodynamically and disrupts the parental molecular 
cloud, producing a classical HII region. Although the evolutionary path of high-mass YSOs 
likely runs through these different categories, this is not yet as well established as for the 
case of low-mass YSOs, and the physical structure and processes dominating each of these 
stages are not well characterized. Herschel observations of water and related species will 
provide an important piece of this puzzle. 



Abundant H2O has been observed with ISO, SWAS and Odin toward the Orion and 
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SKrB2 high-mass st ar-forming; regions (e.g.. ICernicharo et al.ll2006l . iPolehampton et al.ll2007l . 
Persson et al.ll2007l ) . Toward other high-mass YSOs, copious hot gas-phase H^O has been de- 
tecte d in absorption with ISO-SWS (Ivan Dishoeck fc Helmichlll996l . lBoonman fc van Dishoeck 
2003[) but not wit h ISO-LWS. SWAS data indicate generally weak H2O 557 GHz emission 
(jSnell et al.l 120001 ). Analysis of these data for selected high- mass star-forming regions im- 
plies sharp gradients in the water abu ndance, together with strong excitation chang es, but 
their origin is not yet well understood (IBoonman et al.ll2003l . Ivan der Tak et al.ll2006l ). More 
generally, the role of shocks and the relative importance of shocks versus passive heating in 
the inner envelope during the high-mass evolution, either in isolated sources or in clusters, is 
currently unclear. In addition, the relation of these phases to the H2O maser sources needs 
further study. Since individual high-mass sources are at much larger distances (typically 
a few kpc) than their low-mass counterparts and thus spatially unresolved, disentangling 
the various processes requires observations of numerous H2O transitions and high spectral 
resolution, as obtained within WISH using HIFI and FAGS. 

Key questions include: (i) What is the chemistry of the warm gas close to young high- 
mass stars in various evolutionary stages, specifically the distribution of H2O? (ii) What is the 
relative importance of shocks vs. UV for the interaction of the stars with their environment? 
(iii) What are the kinematics of the cold and warm gas close to young high-mass stars? (iv) 
What are the effects of clustered star formation and feedback by protostellar outflows on 
high-mass regions? 

Source and line selection. The WISH Herschel observations of H2O and related species 
cover each of the above phases in the formation of high-mass stars. Pointed observations 
with HIFI are made of a sample of fou r massive pre-stella r core s and 19 pro t ostella r sources. 
The pre-stellar cores are taken from iGarey et al.l (119981 ) and iPillai et al.l ( l2006al ) and are 
observed only in the 557 GHz line down to 10 mK rms. The protostellar list consists of ten 
HMPOs (five mid-infrared bright and five quiet), four hot cores and five U G HII regions. 
The s o urces are nearby (<few kpc) and isolated, selected from surveys by iMolinari et al. 
f ll996h . lSridharan et al.l f l2002l ). IWood fc Ghurchwelll fll989l ). lvan der Tak et al.l J2000bh . The 
basic H2O line list is the same as for the low- mass and intermediate- mass YSOs, but with 
additional Hg^O and Hg^O lines added because the main isotopic lines are very bright and 
highly optically thick. Absorption lines of H2O at high frequencies against the optically thick 
continuum can uniquely probe the kinematics of the protostellar environment (see § 12. 2p . 
Integration times are such that 50-150 mK rms is reached, depending on the line. 

Full PAGS spectral scans of all protostellar sources are taken and provide an unbiased 
overview of low- and high-excitation H2O lines (and thus the shocked ~1000 K gas in the 
beam) together with OH, [O I] and GO at the position of the protostar, as well as spatial 
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information on a 50" scale. S/N ratios >100 on the continuum are reached so that spectrally 
unresolved lines can be detected. To characterize the effects of clustered star formation and 
feedback, including the chemistry and cooling rates of intra-cluster gas, HIFI mapping of 
the P-H2O 1113 GHz line combined with ^^CO 10-9 plus PACS mapping of 0-H2O 1670 
and 1717 GHz down to 300 mK rms is planned over a region of ~ 3' x 3' for six nearby 
cluster-forming clouds where several hot cores and UCHII regions are known to co-exist in 
a small area. 

I nitial results. Figure [H I includes the lowest two p— H2O lines for the high-mass YSO W3 
IRS5 ( Chavarria et al. 201o[ ). The profile of the excited 2o2 — In line reveals the same broad 
and medium-broad components as seen for their low- and intermediate-mass counterparts, 
but the profile of the ground-state line is clearly more complex. In addition to self-absorption 
and absorption of the cold envelope against the continuum, absorptions due to foreground 
clouds lying along the line of sight are seen at various velocities. Moreover, for line-rich 
sources, weak emission lines due to species other than H2O are detected. The ability of HIFI 
to disentangle the different phy sical components is de monstrated by initial performance 
verification data on DR21 (OH) ( Ivan der Tak et al.ll2010l ). Using a slab model and assuming 
an o/p ratio of 3 ( iLis et al.ll2010l ). the broad component due to the outflow is found to have 
an H2O abundance of a few xlO~^. Water is much less abundant in the foreground cloud 
where p— H2O has an abundance of 4 x 10~^. The outer envelope abundance is even lower, 
a few xlO-^°. 

To investigate w hether these abundances are typical of high-mass protostellar sources, 
Marseille et al.l ( I2OIOI ) model four additional high-mass YSOs in the lowest two p— H2O lines 
combined with the ground-state p—B.1^0 hne. Thanks to the high spectral resolution, the 
envelope emissio r i and absorption can be cleanly separated from other physical components. 
Chavarria et al.l ( 120101 ) analyze a larger set of lines for W3 IRS5. Since the same modeling 
philosophy is used for all six sources, the inferred abundances can be directly compared. 
The outer envelope water abundances are low in all cases but show two orders of magnitude 
variation from < 5 x 10~^° to 4 x 10~^. The origin for these variations is not yet understood, 
since no correlation is found with luminosity or evolutionary stage within this small sample. 
For W3 IRS5, a jump in abundance in the inner envelope at T ^ 100 K is inferred from the 
isotopologue lines, the only indication so far of water emission from a hot core. 

Broad emission profiles due to outflows associated with the high-mass protostars are 
commonly seen in the high-mass spectra, testifying to the importance of shocks. Moreover, 
the presence of blue-shifted absorption suggests expansion of the outer envelope rather than 
infall. Finally, the high S/N absorption data reveal cold clouds surrounding the protostellar 
envelope at velocities offset by just a few km s~^, illustrating the complexity of the protostel- 
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lar en vironment and the unique ability of HIFI and water to detect this gas ( iMarseille et aL 



2OIOI). 



In summary, the available data on the first six sources show that questions (i) and (iii) 
can be addressed quantitatively. Questions (ii) and (iv) require more data to answer. 



4.6. Radiation diagnostics 



Motivation. The H2O chemistry can be strongly affected by high-energy radiation from 
YSOs, both UV and X-rays. For typical X-ray luminosities of 10^^ — 10^^ erg s~^ appropriate 
for low-mass YSOs, models show that the H2O abundance decreases by two orders of mag- 
nitude or more in the inn er envelope, where t he X-rays penetrate up to distances of ~1000 
AU (one Herschel beam; IStauber et al.ll2006l ). UV and X-rays also affect the chemistry in 
high- mass YSO envelopes, but to a different degree (IStauber et al.l 120041 : 120051 ). H2O is de- 
stroyed not only by direct photodissociation but also by reactions with ionized molecules, 
so that it is important to understand their chemistry in relation to that of H2O (Fig. |3]). 
Because H2O is not a unique diagnostic of either UV or X-rays, independent diagnostics 
need to be observed. S mall hydrides, especially QH"^, CH, CH+ ari d SH+, are particularly 
powerful in this regard ( IStauber et al.l 120051 . iBruderer et al.l l2010al ). Recent 2D modeling 
has shown that most of the emission from UV tracers originates in the outflow walls, where 
their abundances are enhanced by orders of magnitude due to UV impinging on the walls 
of the outflow cavities l eading to ionization, dissociation and high-temperature chemistry 
jBruderer et"^l2009by : boiOal ). 



Key questions are: (i) Can we establish the presence of UV and/or X-rays through 
chemistry in deeply embedded objects, where the UV or X-rays cannot be observed directly? 
If so, can we quantify it? (ii) How do protostellar UV and X-rays influence the surroundings, 
both the temperature and the water chemistry? 

Source and line selection. Deep pointed observations with HIFI are made of lines of 
seven hydrides that are related to the water network and sensitive to either FUV or X-rays 
(or both) in two high-mass and one low-mass YSO down to 15 mK rms. Subsequently, a 
larger sample of nine low-, intermediate- and high-mass YSOs is surveyed in the lines of CH, 
CH"*" and OH"*". Lines of some hydrides are also obtained serendipitously with the main H2O 
settings in the other subprograms. For example, the CH 536 GHz triplet is often obtained 
together with the o-H2^^0 lio - loi 548 GHz line and the H2O+ In - Oqo 1115 GHz line 
with the p— H2O 111 — O nn 1113 GHz line . Line frequencies have been collected, and in some 
cases first predicted, by iBrudererl (120061 ). 
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Fig. 10. — Collection of lines of hydri des related to th e H2O chemistry detected toward 
the high- mass YSO W3 IRS5 (based on lBenz et al.ll2010l ). The dashed red lines indicate the 
position of the main component at the LSR velocity of the protostellar envelope at —38.4 km 
s~^. For species with fine-structure components (OH"*", H2O"'"), the length of the vertical bars 



indicates the theoretical relative intensities. Note that strong OH"*" and H2O 
occurs around Vlsr ~ km s~^ in diffuse foreground clouds. 
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Initial results. An early surprise from all Herschel-HIFI KPs is the detection of widespread 



of massive YSOs and outflows ( 


Gerin et al. 


2010. 


Ossenkopf et al. 


2010, 


Benz et al. 


2010 




Bruderer et aliboiOb. 


Wvrowski et al.lbolOb 


Neufeld et al. 


2010, 
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2010 


Gupta et al. 



Mkr231 (Ivan der Werf et al.l 120101 ). Model analyses point to an origin in low density gas 
(< 10'^ cm~^) with a low molecular fraction (i.e., a H2/H ratio of a few %) to prevent H2O"'" 
from reacting rapidly with H2 to form HsO^. The low molecular fraction can be the result 
of either high UV or high X-ray fluxes. 

Within WISH, deep integrations on two high-mass YSOs reveal most of the targeted 
molecules, a vindication of the model predictions. CH, NH, OH"*" and H^O"^ are detected 
toward both W3 IRS5 and AFGL 2591 jBenz et allboidlBruderer et al.lboiObh . In addition 
SH+ and H3O+ are seen toward W3 and CH+ toward AFGL 2591. Only SH and NH+ are 
not detected. 1120"*", 0H+ and SH+ are new molecules, with OH"*" and SH"*" o nly rec ently 
detected for the firs t time in ground-based observations by IWyrowski et al.l (l2010al ) and 
Menten et al.l (120111 ). respectively. Figure [TD] contains a collection of spectra of hydrides 
most closely related to water. 

The spectra show a surprising mix of absorption and emission features. Three molecules, 
NH, OH"*" and H2O"'", appear primarily in absorption, either in diffuse foreground clouds along 
the line of sight, in blue-shifted outflows, or in some cases in clouds close to the systemic 
velocity. For these molecules, the J = 1 level lies at 47 K or higher. Limits on the excited 
state lines give rotational temperatures <13 and <19 K for 0H+ and H2O"'", respectively. 
CH and SH"*", for which the J = 1 line lies at less than 26 K, are seen in emission, although 
CH has superposed absorption. The CH+ J=l-0 line occurs mostly in absorption but the 
2-1 line is in emission in AFGL 2591. For W3 IRS5, at least three HsO"*" lines from energy 
levels lying at a few hundred K above ground are detected, which combined with ground- 
based data (jPhillips et al.lll992l ) give a rotational temperature of 240 ± 40 K. Note that all 
excitation temperatures are lower limits to the kinetic temperature of the gas because the 
lines have high critical densities (> 10^ cm~^) resulting in subthermal excitation. 

The observed range in line profiles and excitation conditions immediately indicates that 
the species originate in different parts of the YSO environment. The OH+/H2O"'" abundance 
ratio of > 1 for the blue-shifted gas in both sources indicates low density and high UV fields 
(iGerin et al.l l2010l ). as expected along the cavity walls at large distances from the source. 
The CH, CH "*" and H.-^O"*" ernission hues must originate closer to the protostar, at densities 



> 10^ cm ^ (IBruderer et al.ll2010bl ). Their abundances are roughly consistent with the 2D 



photochemical models, in which high temperature chemistry and high UV fluxes boost their 
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abundances along the outflow walls ( iBruderer et al.ll2009at l2010a] ) . The UV field is enhanced 
by at least four orders of magnitude compared with the general interstellar radiation field in 
these regions. 

Because the H2O+ ground state line occurs close to that of p— H2O, observations of 
both species are available for a l a rge sam ple of high-mass protostars within WISH. Using 
also Hg^O data, IWyrowski et al.l ( l2010bl ) provide column densities for both species for all 
components along 10 lines of sight. H2O"'" is always in absorption, even when H2O is seen 
in emission. Consistent with the above scenario, the highest H2O+ column densities are 
found for the outfiow components. Overall, H2O+/H2O ratios range from 0.01 to > 1, with 
the lower values found in the dense protostellar envelopes and the larger values in diffuse 
foreground clouds and outflows. 

In summary, most of the targeted hydrides are readily detected and they appear to be 
even more abundant and widespread than expected prior to Herschel. When associated with 
protostars, they trace regions of high temperatures and enhanced UV fluxes, such as found 
along cavity outflow walls. However, the fact that some hydrides are seen in absorption 
whereas others appear in emission is not yet fully understood. 



4.7. Circumstellar disks 



Motivation. A fraction of the gas and dust from the protostellar envelope ends up in 
a rotating disk around the pre-main sequence star (< 10 Myr), thereby providing the basic 
material from which planetary systems can be formed. Recent observations and modeling 
have shown that many young disks have a strongly flaring structure in which the surface 
layers are exposed to intense ultraviolet radiation from the star. This results in a highly 
layered structure with a hot top layer which is mostly atomic due to rapid photodissociation; 
an intermediate warm layer with an active chemistry and responsible for most of the observed 
molecular e mission; and a cold, dense midplane lay er where most species are frozen out onto 



grams (e.g.. 



Aikawa et al.ll2002l . lBergin et al.ll2007f ) 



A major question in disk evolution and planet formation studies is the amount of grain 
growth and vertical mixing as a function of radial position and time. In the static mod- 
els described above, httle gas-phase H2O is expected in the outer disk, with most of it 
frozen out as H2O ice in the midplane. However, if the H2O ice is regularly circulated to 
the warmer upper la yers, the gas-phase H2O abundance increases by orders of magnitude 



( iDominik et al. 



2OO5I ). This in turn affects the gas temperature and the abundances of other 



species, including its photodissociation products OH and O. Indeed, analysis of CO observa- 



- 51 - 



tions shows that significa nt gaseous CO is present at temperatures below that for freeze-out 



[e.g., iDartois et al.ll2003[ ). Both photodesorption of ices and rapid vertical mixing combined 
with a steep vertical temperature gradient and grain gro wth have been invoked to explain 
this f lAikawall2007l . loberg et aPboOTi Eersant et al.lboOQl ). 



HIFI observations are uniquely suited to constrain vertical mixing models because H2O 
is the major ice reservoir and its line emission is particularly sensitive to temperature. HIFI 
thus probes the cold H2O gas reservoir in the outer (> 30 AU) part of disks. These ob- 
servations are therefore highly compl ementary to Spitzer detections of hot H2O lines which 
originate in the inner few AU of disks (ISalyk et al.ll2008l . ICarr &: Najitall2008l ). Key questions 
are: (i) What do H2O observations of circumstellar disks tell us about grain growth and the 
importance of vertical mixing from the cold midplane to the warm surface layers where ices 
photodesorb or thermally vaporate? (ii) What is the relative importance of photodissocia- 
tion and desorption? (iii) What are the implications of the Herschel observations of the cold 
water reservoir in disks for the transport of water to the inner planet-forming zones of disks? 

Source and line selection. Deep integrations down to a ~1.5 mK rms in 0.5 km s~^ bin 
are performed for the o— H2O 557 GHz line in four targets (TW Hya, DM Tau, LkCal5 and 
MWC 480) and down to ~4 mK rms for the P-H2O 1113 GHz line in two targets (TW Hya 
and DM Tau). The primary aim of these very deep integrations is to put stringent limits 
on the cold water reservoir. Shorter integrations down to 10 mK rms in the 557 GHz line 
are performed for eight additional targets to determine whether the four deep targets are 
unusual in their H2O emission. The excited H2O 3i2 — 3o3 1097 GHz and CO J=10-9 line 
are observed in two sources to probe the warmer H2O reservoir with HIFI. Sources have 
been selected t o be clean from surrounding cloud emission in single dish observations (cf. 
Thi et all [2001! ). 



Initial results. Deep integrations on the DM Tau disk show a tentative de tection of the 



o— H2O line with a peak temperature of Tmb=2.7 mK and a width of 5.6 km s ^ (IBergin et al. 



2010bl ). Regardless of the reality of this detection and the uncertainty in the coUisional rate 
coefficients, the inferred column density is a factor of 20-130 weaker than predicted by a 
simple model in which gaseous water is produced by photodesorption of icy grains in the UV 
illuminated regions of the disk. The most plausible implication of this lack of water vapor 
is that more than 95-99% of the water ice is locked up in large coagulated grains that have 
settled to the midplane and are not cycled back up to the higher disk layers. Thus, the deep 
HIFI data can directly address the key questions of this subprogram, even without actual 
detections. 
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5. Discussion 

Although less than 10% of the WISH data have been analyzed to date, several general 
trends are emerging. 



5.1. Trends with evolutionary stage 

The development of the spectral features with evolutionary phase is illustrated by Fig. H] 
for low-mass sources. Water emission is not detected from cold dense cores prior to star 
formation, down to very low levels. The implication is that most of the water is frozen out 
as ice on dust grains in the densest regions with at most a small layer of water vapor due to 
photodesorption of ice in the outer regions of the core. As soon as the protostar turns on in 
the center, however, prominent H2O lines appear. One of the main surprises for the low mass 
YSOs is that this emission has broad profiles, even for H2^^0, indicating an origin in shocked 
gas. This association with outflows is directly revealed by imaging of the water line (Fig. |9]). 
H2O abundances in the outflow are high, of order 10~^ — 10~^ (assuming a CO/H2 abundance 
ratio of 10~^), increasing with higher velocity. In contrast, quiescent envelope abundances, 
as derived from narrow absorption features and the lack of narrow H2^^0 emission, are low, 
of order 10~^. 

The last evolutionary stage of protoplanetary disks when the envelope has dissipated 
shows again little H2O emission down to very low levels (Fig. Hj). This is the first time that 
the cold water reservoir in the outer regions of disks has been probed. The lack of emission 
implies that the upper layers of this disk are surprisingly 'dry', both in the gas and in the 
solid phase. About 95-99% of the water ice appears to be contained in large icy grains in the 
disk midplane which are not circled back up to higher layers. If confirmed by observations 
of additional sources, this may have profound implications for our understanding of the 
chemistry and physics of the outer disks. On th e other hand, lines of war m and hot water 



have been detected in other disks with Spitzer (jPontoppidan et al.l I2OIOI ) and with PACS 



(ISturm et al.ll2010l ). although the latter detection needs confirmation.. Tying these data 
together should provide a nearly complete view of the distribution of water vapor from the 
inner to the outer disk. 

For all of these stages, the availability of data on both o— H2O and p— H2O indicates 
that the low water abundance is robust and independent of any o/p ratio of water. The 
overall emerging picture is therefore that most water is formed on grains in the molecular 
cloud and pre-stellar core phases, is present mostly as ice in the quiescent envelope and is 
likely transported in this form into the outer part of young disks without further evaporation. 
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consistent with evolutionary models by IVisser et al.l (120091 ) 



5.2. Trends with luminosity 

Fig. [6] illustrates the trend with luminosity for the embedded protostellar evolution 
stage. An initial conclusion is that the line profiles are remarkably similar in spite of the 
four orders of magnitude range in luminosity. All profiles show a broad (FWHM 25-40 km 
s~^) component due to the shocks along the outflow cavity walls on scales > 1000 AU ('shell 
shocks'), and a medium-broad (FWHM 5-10 km s~^) component ascribed to shocks caused 
by the interaction of jets and winds with the inner dense (> 10^ cm~^) envelope on scales 
< 1000 AU. Thus, the processes by which the jets and winds that are launched close to the 
star interact with their environment are apparently similar across a wide range of masses 
and luminosities. The H2O medium-broad profiles clearly have larger line widths than those 
of other high density envelope tracers such as HCO"*" 4-3 and CS 7-6, but generally agree 
with the profiles of other grain surface products such as CH3OH. 

For low- and intermediate mass YSOs, the outflows dominate the H2O emission whereas 
for high-mass YSOs, narrower lines from the passively heated envelope become more promi- 
nent, especially for 132^*^0 and H2^^0. Inferred outer envelope abundances based on the 
isotopologue data and narrow absorption features vary between 10~^ and lO"^*' with no ob- 
vious trends across luminosity or evolutionary stage, except for a potentially higher value 
for one intermediate- mass source. These low abundances thus confirm the earlier SWAS 
conclusions but can significantly refine them because of the availability of more lines from 
both o— and p— H2O allowing more accurate determinations. Because of the dominance 
of outflow emission, it is more difficult than expected to put firm constraints on the inner 
quiescent 'hot core' water abundances. Generally the 'stirred-up' inner region probed by 
the medium-broad component overwhelms any hot core emission, even though this warm 
quiescent gas is clearly seen in, for example, C^^O J=9-8 or 10-9 emission. Simultaneous 
analysis of these CO high— J data with the H2O lines may help to further disentangle the 
components. Indeed, the coverage of the entire CO ladder up to J = 35 with all three instru- 
ments on Herschel will soon become one of the most important diagnostics to characterize 
the physical processes in a wide variety of sources. 

The PACS and HIFI observations of high-J CO, H2O and [O I] illustrate the need for 
the development of multi-dimensional models which include UV photon and shock heating 
along the cavity walls for both low- and high-mass sources. This need is strengthened by the 
detection of various hydrides enhanced by UV radiation, including OH"*", H2O"'", SH"*", CH 
and CH"*". Initial quantitative estimates of the total far-infrared line emission indicate that 
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H2O is not the major coolant of the gas, even in the densest (> 10^ cm~^) regions, but this 
needs to be confirmed by observational data at more positions, both on and off source. [O 
I] and OH have higher luminosities and appear to trace a dissociative shock in the densest 
inner part. 

5.3. Water chemistry 

The Herschel data demonstrate that all three chemical networks illustrated in Fig. [3] are 
important in star-forming regions. Solid-state chemistry dominates in the coldest regions, 
as evidenced by the absence of H2O line emission from cold pre-stellar cores and the low 
gaseous water abundances in the outer protostellar envelopes. Most of the water is in the 
form of ice formed through the reactions illustrated in the bottom part of Fig. [31 The 
tiny amounts of water gas detected with HIFI either result from non-thermal desorption 
of water ice or from gas-phase ion-molecule reactions, or from a combination of both. If 
the absorption of H2O gas toward the pre-stellar core L1544 is confirmed, this first direct 
detection of water vapor in a quiescent dark cloud will allow a quantitative test of these 
chemical models. Similarly, a thorough analysis of the abundance structure of H2O gas and 
ice for those protostellar envelopes where both species are detected can test key aspects of 
the cold chemistry. Eventually, fully coupled time-dependent gas-grain chemistry models 
should include the new solid-state laboratory data on water ice formation. Such models can 
then also be used to interpret the puzzling HDO/H2O abundance ratios found in different 
regions. 

The detection of all the hydrides involved in the ion-molecule scheme leading to water 
— OH"*", H2O+, H3O+, HCO"*" and OH — will be one of the lasting legacies of the Herschel 
mission. While these species have been part of the reaction networks for decades, the ease 
with which rapidly reacting ions such as 0H+ and H2O"'" are detected has surprised even 
the astrochemists. Since all key reaction rates involved are well understood, the abundance 
ratios can be used to probe the physics of the regions in which they are found, such as high 
UV fields and low densities. The final clue to the basic oxygen gas-phase chemistry network 
will be provided by deep searches for O2, which are planned for some of the WISH sources 
in other key programs. 

The bright water emission from protostars is dominated by C— type shocks where the 
temperature rises to a few thousand K. The high-quality line profiles may allow the two 
chemical regimes to be separated: sputtering of icy mantles and high-temperature chemistry. 
Whether all available oxygen is driven into water at the highest temperatures and velocities is 
still an open question. OH/H2O ratios derived from velocity-resolved OH line profiles could 
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form a powerful test of the neutral-neutral reactions. If the contribution from dissociative 

J— type shocks can be cleanly separated, as appears to be the case for HH 46, the observed 
0/OH abundance ratios would provide a further test of the high temperature gas-phase 
reaction networks, but now under conditions with much higher H/H2 ratios in the presence 
of UV radiation. 

Finally, all three chemical networks play a role in the water chemistry in disks, with 
high temperature chemistry dominating in the inner disk and ice chemistry in the outer disk. 
Here Herschel teaches us another lesson: the importance of coupling the water chemistry 
with the dynamics and physical evolution of the disks, in particular the growth and settling 
of the icy grains to the midplane. 

6. Conclusions 

The goal of the WISH program is to use water and related species as physical and 
chemical probes of star-forming regions over a range of luminosities and evolutionary stages. 
The initial results presented in §4 and discussed in §5 demonstrate that many of the questions 
asked in the individual subprograms can indeed be addressed by the WISH data, thanks to 
the combination of excellent sensitivity, fully resolved line profiles and spatial information, 
thereby validating the observational strategy. Initial surprises include a near absence of 
gaseous water in pre-stellar cores and disks, the dominance of shocks rather than hot cores 
in controlling the bright water emission from protostars, and the detection of all of the ions 
and hydrides involved in the water chemistry schemes. For pre-stellar cores our observations 
are in agreement with models that suggest near total freeze-out of water in the form of ice, 
whereas in disks the inferred water abundance is at levels significantly lower than predicted. 
So far, water has not yet been found to be the major coolant in any of these regions. All 
remaining observations are expected to be taken within the next year. Quantitative analysis 
will require the further development of multidimensional models of protostellar envelopes, 
outflows and disks. Together with the results from related Herschel key programs, they will 
greatly enhance our understanding of water in the galactic interstellar medium and solar 
system, and provide a true legacy to follow the water trail from the most diffuse gas to dense 
cores and disks, and eventually comets and planets in our own Solar system. 



-56- 



7. Acknowledgments 

The authors are grateful to the HIFI project scientist, Xander Tielens, for helping to 
ensure the powerful HIFI capabilities for water observations, and for his encouragements 
over the last decade. They salute the HIFI and PACS instrument builders for providing 
two superb scientific instruments, and are much indebted to the laboratory and theoretical 
chemistry groups for providing the necessary molecular data to analyze and interpret the 
water data. They thank the referee, David HoUenbach, for his constructive comments on the 
manuscript and Laurent Wiesenfeld for helpful discussions on the H2O-H2 collisional rate 
coefficients. We thank the Netherlands Organization for Scientific Research (NWO) and 
many other national funding agencies for their financial support. 

HIFI has been designed and built by a consortium of institutes and university de- 
partments from across Europe, Canada and the United States under the leadership of 
SRON Netherlands Institute for Space Research, Groningen, The Netherlands and with ma- 
jor contributions from Germany, France and the US. Consortium members are: Canada: 
CSA, U.Waterloo; France: CESR, LAB, LERMA, IRAM; Germany: KOSMA, MPIfR, 
MPS; Ireland, NUI Maynooth; Italy: ASI, IFSLINAF, Osservatorio Astrofisico di Arcetri- 
INAF; Netherlands: SRON, TUD; Poland: CAMK, CBK; Spain: Observatorio Astronomico 
Nacional (IGN), Centro de Astrobiologia (CSIC-INTA). Sweden: Chalmers University of 
Technology-MC2, RSS & GARD; Onsala Space Observatory; Swedish National Space Board, 
Stockholm University - Stockholm Observatory; Switzerland: ETH Zurich, FHNW; USA: 
Caltech, JPL, NHSC. 

HCSS / HSpot / HIPE are joint developments by the Herschel Science Ground Segment 
Consortium, consisting of ESA, the NASA Herschel Science Center, and the HIFI, PACS 
and SPIRE consortia. 



-57- 



REFERENCES 

Aikawa, Y. 2007, ApJ, 656, L93 

Aikawa, Y., Herbst, E., Roberts, H., & Caselli, P. 2005, ApJ, 620, 330 

Aikawa, Y., van Zadelhoff, G. J., van Dishoeck, E. F., & Herbst, E. 2002, A&A, 386, 622 

Aikawa, Y., Wakelam, V., Garrod, R. T., & Herbst, E. 2008, ApJ, 674, 984 

Andre, P., & Montmerle, T. 1994, ApJ, 420, 837 

Andre, P., Ward-Thompson, D., & Barsony, M. 2000, in Protostars and Planets IV, ed. 
V. Mannings, A. P. Boss, S.S. Russell (Tucson: Univ. of Arizona), 59 

Atkinson, R., Baulch, D., Cox, R., Crowley, J., Hampson, R., & Hynes, R. 2004, Atmospheric 
chemistry and physics, 4, 1461 

Bachiller, R., Codclla, C, Colomer, F., Liechti, S., & Walmsley, C. M. 1998, A&A, 335, 266 

Bachiller, R., & Tafalla, M. 1999, in NATO ASIC Proc. 540: The Origin of Stars and 
Planetary Systems, ed. C. J. Lada & N. D. Kylafis, 227 

Beltran, M. T., Cesaroni, R., Neri, R., Codella, C, Furuya, R. S., Testi, L., & Olmi, L. 2005, 
A&A, 435, 901 

Benz, A. O., Stauber, P., Bourke, T. L., van der Tak, F. F. S., van Dishoeck, E. F., & 
J0rgensen, J. K. 2007, A&A, 475, 549 

Benz, A. O., et al. 2010, A&A, 521, L35 

Bergin, E. A., Aikawa, Y., Blake, G. A., & van Dishoeck, E. F. 2007, in Protostars and 
Planets V, ed. B. Reipurth, D. Jewitt, K. Keil (Tucson: Univ. of Arizona), 751-766 

Bergin, E. A., Alves, J., Huard, T., & Lada, C. J. 2002, ApJ, 570, LlOl 

Bergin, E. A., & Snell, R. L. 2002, ApJ, 581, L105 

Bergin, E. A., & Tafalla, M. 2007, ARA&A, 45, 339 

Bergin, E. A., et al. 2000, ApJ, 539, L129 

— . 2010a, A&A, 521, L20 

— . 2010b, A&A, 521, L33 



-58- 

Beuther, H., Schilke, P., Gueth, F., McCaughrean, M., Andersen, M., Sridharan, T. K., & 
Menten, K. M. 2002a, A&A, 387, 931 

Beuther, H., Schilke, P., Menten, K. M., Motte, F., Sridharan, T. K., & Wyrowski, F. 2002b, 
ApJ, 566, 945 

Beuther, H., Walsh, A. J., Thorwirth, S., Zhang, Q., Hunter, T. R., Megeath, S. T., & 
Menten, K. M. 2008, A&A, 481, 169 

Beuther, H., Zhang, Q., Bergin, E. A., Sridharan, T. K., Hunter, T. R., & Leurini, S. 2007, 
A&A, 468, 1045 

Bjerkeh, P., et al. 2009, A&A, 507, 1455 

Bockelee-Morvan, D., et al. 1998, Icarus, 133, 147 

Boogert, A. C. A., et al. 2008, ApJ, 678, 985 

Boonman, A. M. S., Doty, S. D., van Dishoeck, E. F., Bergin, E. A., Melnick, G. J., Wright, 
C. M., & Stark, R. 2003, A&A, 406, 937 

Boonman, A. M. S., & van Dishoeck, E. F. 2003, A&A, 403, 1003 

Brinch, C., Crapsi, A., Hogerheijde, M. R., & J0rgensen, J. K. 2007, A&A, 461, 1037 

Brinch, C., & Hogerheijde, M. R. 2010, A&A, 523, A25 

Bronfman, L., Nyman, L., & May, J. 1996, A&AS, 115, 81 

Bruderer, S. 2006, Master's thesis. Institute of Astronomy, ETH Zurich, Switzerland 

Bruderer, S., Benz, A. O., Doty, S. D., van Dishoeck, E. F., & Bourke, T. L. 2009a, ApJ, 
700, 872 

Bruderer, S., Benz, A. O., Stauber, P., & Doty, S. D. 2010a, ApJ, 720, 1432 
Bruderer, S., Doty, S. D., & Benz, A. O. 2009b, ApJS, 183, 179 
Bruderer, S., et al. 2010b, A&A, 521, L44 

Butner, H. M., Evans, II, N. J., Harvey, P. M., Mundy, L. G., Natta, A., & Randich, M. S. 
1990, ApJ, 364, 164 

Butner, H. M., Evans, II, N. J., Lester, D. F., Levreault, R. M., & Strom, S. E. 1991, ApJ, 
376, 636 



-59- 

Carey, S. J., Clark, F. O., Egan, M. P., Price, S. D., Shipman, R. F., & Kuchar, T. A. 1998, 
ApJ, 508, 721 

Carr, J. S., & Najita, J. R. 2008, Science, 319, 1504 

Caselli, P., Walmsley, C. M., Tafalla, M., Dore, L., & Myers, P. C. 1999, ApJ, 523, L165 

Caselli, P., Walmsley, C. M., Zucconi, A., Tafalla, M., Dore, L., & Myers, P. C. 2002, ApJ, 
565, 331 

Caselli, P., et al. 2010, A&A, 521, L29 

Ceccarelli, C, Castets, A., Caux, E., HoUenbach, D., Loinard, L., Molinari, S., & Tielens, 
A. G. G. M. 2000a, A&A, 355, 1129 

Ceccarelli, C, HoUenbach, D. J., & Tielens, A. G. G. M. 1996, ApJ, 471, 400 

Ceccarelli, C, Loinard, L., Castets, A., Tielens, A. G. G. M., & Caux, E. 2000b, A&A, 357, 
L9 

Ceccarelli, C, et al. 1999, A&A, 342, L21 
— . 2010, A&A, 521, L22 

Cernicharo, J., & Crovisier, J. 2005, Space Science Reviews, 119, 29 

Cernicharo, J., Thum, C, Hein, H., John, D., Garcia, P., & Mattioco, F. 1990, A&A, 231, 
L15 

Cernicharo, J., et al. 2006, ApJ, 649, L33 

Cesaroni, R., Hofner, P., Walmsley, C. M., & Churchwell, E. 1998, A&A, 331, 709 

Cesaroni, R., Neri, R., Olmi, L., Testi, L., Walmsley, C. M., & Hofner, P. 2005, A&A, 434, 
1039 

Chandler, C. J., Brogan, C. L., Shirley, Y. L., & Loinard, L. 2005, ApJ, 632, 371 
Charnley, S. B. 1997, ApJ, 481, 396 
Chavarria, L., et al. 2010, A&A, 521, L37 

Chen, J., Evans, N. J., Lee, J., & Bourke, T. L. 2009, ApJ, 705, 1160 

Chen, X., Launhardt, R., Bourke, T. L., Henning, T., & Barnes, P. J. 2008, ApJ, 683, 862 



Cheung, A. C, Rank, D. M., & Townes, C. H. 1969, Nature, 221, 626 
Codella, C, et al. 2010, A&A, 518, L112 

Crapsi, A., Caselli, P., Walmsley, C. M., Myers, P. C, Tafalla, M., Lee, C. W., & Bourke, 
T. L. 2005, ApJ, 619, 379 

Crapsi, A., van Dishoeck, E. F., Hogerheijde, M. R., Pontoppidan, K. M., & DuUemond, 
C. P. 2008, A&A, 486, 245 

Crimier, N., Ceccarelli, C, Maret, S., Bottinelli, S., Caux, E., Kahane, C, Lis, D. C, & 
Olofsson, J. 2010, A&A, 519, A65 

Cuppen, H. M., loppolo, S., Romanzin, C, & Linnartz, H. 2010, Phys. Chem. Chem. Phys., 
12, 12077 

Dalgarno, A., & Black, J. H. 1976, Reports on Progress in Physics, 39, 573 

Daniel, P., Dubernet, M., Pacaud, P., & Grosjean, A. 2010, A&A, 517, A13 

Dartois, E., Thi, W., Geballe, T. R., Deboffle, D., d'Hendecourt, L., & van Dishoeck, E. 
2003, A&A, 399, 1009 

De Buizer, J. M., & Minier, V. 2005, ApJ, 628, L151 

de Geus, E. J., Bronfman, L., & Thaddeus, P. 1990, A&A, 231, 137 

de Graauw, T., et al. 2010, A&A, 518, L6 

Di Francesco, J., Johnstone, D., Kirk, H., MacKenzie, T., & Ledwosinska, E. 2008, ApJS, 
175, 277 

Dick, M., Drouin, B., & Pearson, J. 2010, Physical review. A, Atomic, molecular, and optical 
physics, 81 

Dominik, C, CeccareUi, C, HoUenbach, D., & Kaufman, M. 2005, ApJ, 635, L85 
Doty, S. D., & Neufeld, D. A. 1997, ApJ, 489, 122 

Doty, S. D., Schoier, F. L., & van Dishoeck, E. F. 2004, A&A, 418, 1021 
Doty, S. D., van Dishoeck, E. F., & Tan, J. C. 2006, A&A, 454, L5 

Doty, S. D., van Dishoeck, E. F., van der Tak, F. F. S., & Boonman, A. M. S. 2002, A&A, 
389, 446 



-61 - 

Dubernet, M., Daniel, F., Grosjean, A., & Lin, C. Y. 2009, A&A, 497, 911 
Dubernet, M., et al. 2006, A&A, 460, 323 

Dulieu, F., Amiaud, L., Congiu, E., Fillion, J., Matar, E., Momeni, A., Pirronello, V., & 
Lemaire, J. L. 2010, A&A, 512, A30 

DuUemond, C. P., & Dominik, C. 2004, A&A, 417, 159 

Dzib, S., Loinard, L., Mioduszewski, A. J., Boden, A. F., Rodriguez, L. F., & Torres, R. M. 
2010, ApJ, 718, 610 

Eiroa, C, Palacios, J., & Casali, M. M. 1998, A&A, 335, 243 

Elitzur, M., & Watson, W. D. 1978, ApJ, 222, L141 

Evans, N. J., et al. 2009, ApJS, 181, 321 

Evans, II, N. J., et al. 2003, PASP, 115, 965 

Fich, M., et al. 2010, A&A, 518, L86 

Franklin, J., Snell, R. L., Kaufman, M. J., Melnick, G. J., Neufeld, D. A., HoUenbach, D. J., 
& Bergin, E. A. 2008, ApJ, 674, 1015 

Eraser, H. J., CoUings, M. P., McCoustra, M. R. S., & Williams, D. A. 2001, MNRAS, 327, 
1165 

Froebrich, D. 2005, ApJS, 156, 169 

Fuente, A., Neri, R., & Caselli, P. 2005a, A&A, 444, 481 

Fuente, A., Rizzo, J. R., Caselli, P., Bachiller, R., & Henkel, C. 2005b, A&A, 433, 535 
Garay, G., Brooks, K. J., Mardones, D., Norris, R. P., & Burton, M. G. 2002, ApJ, 579, 678 
Gensheimer, P. D., Mauersberger, R., & Wilson, T. L. 1996, A&A, 314, 281 
Gerin, M., et al. 2010, A&A, 518, LllO 

Giannini, T., McCoey, C, Nisini, B., Cabrit, S., Caratti o Garatti, A., Calzoletti, L., & 
Flower, D. R. 2006, A&A, 459, 821 

Giannini, T., Nisini, B., & Lorenzetti, D. 2001, ApJ, 555, 40 

Giannini, T., et al. 2005, A&A, 433, 941 



-62 - 

Gibb, A. G., Wyrowski, F., & Mundy, L. G. 2004a, ApJ, 616, 301 

Gibb, E. L., Whittet, D. C. B., Boogert, A. C. A., & Tielens, A. G. G. M. 2004b, ApJS, 
151, 35 

Goldsmith, P. F., & Langer, W. D. 1978, ApJ, 222, 881 
Green, S., Maluendes, S., & McLean, A. D. 1993, ApJS, 85, 181 

GuUixson, C, Gehrz, R. D., Hackwell, J. A., Grasdalen, G. L., & Castelaz, M. 1983, ApJS, 
53, 413 

Gupta, H., et al. 2010, A&A, 521, L47 

Gustafsson, M., Ravkilde, T., Kristensen, L. E., Cabrit, S., Field, D., & Pineau Des Forets, 
G. 2010, A&A, 513, A5 

Hachisuka, K., et al. 2006, ApJ, 645, 337 

Hartogh, P., et al. 2010, A&A, 518, L150 

Harwit, M., Neufeld, D. A., Melnick, G. J., & Kaufman, M. J. 1998, ApJ, 497, L105 
Helmich, F. P., et al. 1996, A&A, 315, L173 
Herbst, E., & Klemperer, W. 1973, ApJ, 185, 505 
Herbst, E., & van Dishoeck, E. F. 2009, ARA&A, 47, 427 

Hersant, F., Wakelam, V., Dutrey, A., Guilloteau, S., & Herbst, E. 2009, A&A, 493, L49 
Hirota, T., et al. 2008, PAS J, 60, 37 
Hjalmarson, A., et al. 2003, A&A, 402, L39 

Hogerheijde, M. R., & van der Tak, F. F. S. 2000, A&A, 362, 697 

Hogerlieijdc, M. R., van Dishoeck, E. F., Blake, G. A., & van Langevelde, H. J. 1997, ApJ, 
489, 293 

Hogerheijde, M. R., van Dishoeck, E. F., Salverda, J. M., & Blake, G. A. 1999, ApJ, 513, 
350 

HoUenbach, D., Kaufman, M. J., Bergin, E. A., & Melnick, G. J. 2009, ApJ, 690, 1497 
HoUenbach, D., & McKee, C. F. 1989, ApJ, 342, 306 



-63- 

Hiibers, H. W., Evenson, K. M., Hill, C, & Brown, J. M. 2009, J. Chem. Phys., 131, 034311 

Hunter, T. R., Brogan, C. L., Indebetouw, R., & Cyganowski, C. J. 2008, ApJ, 680, 1271 

loppolo, S., Cuppen, H. M., Romanzin, C, van Dishoeck, E. P., & Linnartz, H. 2008, ApJ, 
686, 1474 

loppolo, S., Cuppen, H. M., Romanzin, C, van Dishoeck, E. P., & Linnartz, H. 2010, Phys. 
Chem. Chem. Phys, 12, 12065 

Ivezic, Z., & Elitzur, M. 1997, MNRAS, 287, 799 

Jacq, T., Walmsley, C. M., Henkel, C, Baudry, A., Mauersberger, R., & Jewell, P. R. 1990, 
A&A, 228, 447 

Johnstone, D., Fich, M., Mitchell, G. P., & Moriarty-Schieven, G. 2001, ApJ, 559, 307 
Johnstone, D., et al. 2010, A&A, 521, L41 
J0rgensen, J. K. 2004, A&A, 424, 589 

J0rgensen, J. K., Bourke, T. L., Myers, P. C, Schoier, P. L., van Dishoeck, E. P., & Wilner, 
D. J. 2005a, ApJ, 632, 973 

J0rgensen, J. K., Johnstone, D., van Dishoeck, E. P., & Doty, S. D. 2006, A&A, 449, 609 

J0rgensen, J. K., Schoier, P. L., & van Dishoeck, E. P. 2002, A&A, 389, 908 

— . 2004, A&A, 416, 603 

— . 2005b, A&A, 437, 501 

J0rgensen, J. K., & van Dishoeck, E. P. 2010a, ApJ, 725, L172 
— . 2010b, ApJ, 710, L72 

J0rgenscn, J. K., van Dishoeck, E. P., Visser, R., Bourke, T. L., Wilner, D. J., Lommen, D., 
Hogerheijde, M. R., & Myers, P. C. 2009, A&A, 507, 861 

J0rgensen, J. K., et al. 2007, ApJ, 659, 479 

Kaufman, M. J., & Neufeld, D. A. 1996, ApJ, 456, 611 

Keene, J., Philhps, T. G., & van Dishoeck, E. P. 1997, in lAU Symposium, Vol. 170, , 
382-384 



Keto, E., & Caselli, P. 2010, MNRAS, 402, 1625 

Keto, E., Rybicki, G. B., Bergin, E. A., & Plume, R. 2004, ApJ, 613, 355 

Klippenstein, S. J., Georgievskii, Y., & McCall, B. J. . 2010, J. Phys. Chem. A, 114, 278 

Knude, J., & Hog, E. 1998, A&A, 338, 897 

Kristensen, L. E., van Dishocck, E. P., van Kempen, T. A., Cuppen, H. M., Brinch, C, 
J0rgensen, J. K., & Hogerheijde, M. R. 2010a, A&A, 516, A57 

Kristensen, L. E., et al. 2010b, A&A, 521, L30 

Ladd, E. P., Deane, J. R., Sanders, D. B., & Wynn- Williams, C. G. 1993, ApJ, 419, 186 
Lee, J., Bergin, E. A., & Evans, II, N. J. 2004, ApJ, 617, 360 
Lefioch, B., et al. 2010, A&A, 518, L113 
Lis, D. C, et al. 2010, A&A, 521, L26 

Liseau, R., Lorenzetti, D., Nisini, B., Spinoglio, L., & Moneti, A. 1992, A&A, 265, 577 
Liu, P., & et al. 2010, A&A in press 

Loinard, L., Castets, A., Ceccarelli, C, Caux, E., & Tielens, A. G. G. M. 2001, ApJ, 552, 
L163 

Lommen, D., J0rgensen, J. K., van Dishoeck, E. P., & Crapsi, A. 2008, A&A, 481, 141 

Maret, S., Ceccarelli, C., Caux, E., Tielens, A. G. G. M., & Castets, A. 2002, A&A, 395, 
573 

Maret, S., Ceccarelli, C, Tielens, A. G. G. M., Caux, E., Lefloch, B., Paure, A., Castets, A., 
& Plower, D. R. 2005, A&A, 442, 527 

Marseille, M., Bontemps, S., Herpin, P., van der Tak, P. P. S., & Purcell, C. R. 2008, A&A, 
488, 579 

Marseille, M. G., et al. 2010, A&A, 521, L32 

Meeus, G., Waters, L. B. P. M., Bouwman, J., van den Ancker, M. E., Waelkens, C, & 
Malfait, K. 2001, A&A, 365, 476 

Melnick, G. J. 2009, in Astronomical Society of the Pacific Conference Series, Vol. 417, , 59 



-65 - 

Melnick, G. J., Tolls, V., Neufeld, D. A., Yuan, Y., Sonnentrucker, P., Watson, D. M., 
Bergin, E. A., & Kaufman, M. J. 2008, ApJ, 683, 876 

Melnick, G. J., et al. 2000, ApJ, 539, L77 

Menten, K. M., Wyrowski, F., Belloche, A., Giisten, R., Dedes, L., & Miiller, H. S. P. 2011, 
A&A, 525, A77 

Miyauchi, N., Hidaka, H., Chigai, T., Nagaoka, A., Watanabe, N., & Kouchi, A. 2008, Chem. 
Phys. Lett., 456, 27 

Molinari, S., Brand, J., Cesaroni, R., & Palla, F. 1996, A&A, 308, 573 

Mookerjea, B., Casper, E., Mundy, L. G., & Looney, L. W. 2007, ApJ, 659, 447 

Motte, F., Bontemps, S., Schilke, P., Schneider, N., Menten, K. M., & Broguiere, D. 2007, 
A&A, 476, 1243 

Motte, F., Schilke, P., & Lis, D. C. 2003, ApJ, 582, 277 

Mozurkewich, D., Schwartz, P. R., & Smith, H. A. 1986, ApJ, 311, 371 

Mueller, K. E., Shirley, Y. L., Evans, 11, N. J., & Jacobson, H. R. 2002, ApJS, 143, 469 

Miiller, H. S. P., Thorwirth, S., Roth, D. A., & Winnewisser, G. 2001, A&A, 370, L49 

Murtz, P., Zink, L. R., Evenson, K. M., & Brown, J. M. 1998, J. Chem. Phys., 109, 9744 

Neufeld, D. A., & Dalgarno, A. 1989, ApJ, 344, 251 

Neufeld, D. A., & Kaufman, M. J. 1993, ApJ, 418, 263 

Neufeld, D. A., et al. 2009, ApJ, 706, 170 

— . 2010, A&A, 521, LIO 

Nisini, B., Antoniucci, S., Giannini, T., & Lorenzetti, D. 2005, A&A, 429, 543 

Nisini, B., Benedettini, M., Giannini, T., Codella, C, Lorenzetti, D., di Giorgio, A. M., & 
Richer, J. S. 2000, A&A, 360, 297 

Nisini, B., Codella, C, Giannini, T., Santiago Garcia, J., Richer, J. S., Bachiller, R., & 
Tafalla, M. 2007, A&A, 462, 163 

Nisini, B., Giannini, T., & Lorenzetti, D. 2002, ApJ, 574, 246 



Nisini, B., et al. 2010, A&A, 518, L120 
Noriega-Crespo, A., et al. 2004, ApJS, 154, 352 

Oberg, K. I., Fuchs, G. W., Awad, Z., Fraser, H. J., Schlemmer, S., van Dishoeck, E. F., & 
Linnartz, H. 2007, ApJ, 662, L23 

Oberg, K. I., Linnartz, H., Visser, R., & van Dishoeck, E. F. 2009, ApJ, 693, 1209 

Olberg, M. 2010, HIFI-ICC Technical note version 1.1, 2010-11-17 

Olmi, L., Cesaroni, R., & Walmsley, C. M. 1996, A&A, 307, 599 

Ossenkopf, V., et al. 2010, A&A, 518, LUl 

Parise, B., Castets, A., Hcrbst, E., Caux, E., Ceccarelli, C, Mukhopadhyay, I., & Tielens, 
A. G. G. M. 2004, A&A, 416, 159 

Parise, B., CeccareUi, C., & Maret, S. 2005a, A&A, 441, 171 

Parise, B., Simon, T., Caux, E., Dartois, E., Ceccarelli, C, Rayner, J., & Tielens, A. G. G. M. 
2003, A&A, 410, 897 

Parise, B., et al. 2005b, A&A, 431, 547 

Persson, C. M., et al. 2007, A&A, 476, 807 

Phillips, T. G., van Dishoeck, E. F., & Keene, J. 1992, ApJ, 399, 533 
Phillips, T. R., Maluendes, S., & Green, S. 1996, ApJS, 107, 467 

Pickett, H. M., Poynter, I. R. L., Cohen, E. A., Delitsky, M. L., Pearson, J. C, & MuUer, 
H. S. P. 1998, J. Quant. Spec. Radiat. Transf., 60, 883 

Pilbratt, G. L., et al. 2010, A&A, 518, LI 

Pillai, T., Wyrowski, F., Carey, S. J., & Menten, K. M. 2006a, A&A, 450, 569 

Pillai, T., Wyrowski, F., Menten, K. M., & Kriigel, E. 2006b, A&A, 447, 929 

Plume, R., et al. 2004, ApJ, 605, 247 

Poelman, D. R., & Spaans, M. 2005, A&A, 440, 559 

Poelman, D. R., & van der Tak, F. F. S. 2007, A&A, 475, 949 



Poglitsch, A., et al. 2010, A&A, 518, L2 

Polehampton, E. T., Baluteau, J., Swinyard, B. M., Goicoechea, J. R., Brown, J. M., White, 
G. J., Cernicharo, J., & Grundy, T. W. 2007, MNRAS, 377, 1122 

Pontoppidan, K. M., Salyk, C, Blake, G. A., Meijerink, R., Carr, J. S., & Najita, J. 2010, 
ApJ, 720, 887 

Pontoppidan, K. M., van Dishoeck, E. P., & Dartois, E. 2004, A&A, 426, 925 
Raymond, S. N., Quinn, T., & Lunine, J. I. 2004, Icarus, 168, 1 
Roberts, H., Herbst, E., & Millar, T. J. 2003, ApJ, 591, L41 
Roberts, H., & Millar, T. J. 2007, A&A, 471, 849 

Robitaille, T. P., Whitney, B. A., Indebetouw, R., Wood, K., & Denzmore, P. 2006, ApJS, 
167, 256 

Rodgers, S. D., & Charnley, S. B. 2003, ApJ, 585, 355 

Rodon, J. A., Beuther, H., Megeath, S. T., & van der Tak, F. F. S. 2008, A&A, 490, 213 
Romanzin, C., , & et al. 2010, MNRAS in press 

Sakai, N., Sakai, T., Hirota, T., Burton, M., & Yamamoto, S. 2009, ApJ, 697, 769 

Salyk, C., Pontoppidan, K. M., Blake, G. A., Lahuis, F., van Dishoeck, E. F., & Evans, II, 
N. J. 2008, ApJ, 676, L49 

Sandell, G. 2000, A&A, 358, 242 

Sandell, G., & Sievers, A. 2004, ApJ, 600, 269 

Schilke, P., et al. 2010, A&A, 521, Lll 

Schoier, F. L., J0rgensen, J. K., van Dishoeck, E. F., & Blake, G. A. 2002, A&A, 390, 1001 
— . 2004, A&A, 418, 185 

Schoier, F. L., van der Tak, F. F. S., van Dishoeck, E. F., & Black, J. H. 2005, A&A, 432, 
369 

Schulz, A., Guesten, R., Walmsley, C. M., & Serabyn, E. 1991, A&A, 246, L55 
Seidensticker, K. J., & Schmidt-Kaler, T. 1989, A&A, 225, 192 



Shevchenko, V. S., & Yakubov, S. D. 1989, Soviet Ast., 33, 370 
Shu, F. H., Adams, F. C, & Lizano, S. 1987, ARA&A, 25, 23 
Slawson, R. W., & Reed, B. C. 1988, AJ, 96, 988 

Snell, R. L., HoUenbach, D., Howe, J. E., Neufeld, D. A., Kaufman, M. J., Melnick, G. J., 
Bergin, E. A., & Wang, Z. 2005, ApJ, 620, 758 

Snell, R. L., Scoville, N. Z., Sanders, D. B., & Erickson, N. R. 1984, ApJ, 284, 176 

Snell, R. L., et al. 2000, ApJ, 539, LlOl 

SoUins, R K., et al. 2004, ApJ, 616, L35 

Sridharan, T. K., Beuther, H., Schilke, R, Menten, K. M., & Wyrowski, F. 2002, ApJ, 566, 
931 

Stanke, T., McCauglirean, M. J., & Zinnecker, H. 2000, A&A, 355, 639 

Stauber, P., Doty, S. D., van Dishoeck, E. F., & Benz, A. O. 2005, A&A, 440, 949 

Stauber, P., Doty, S. D., van Dishoeck, E. F., J0rgensen, J. K., & Benz, A. O. 2004, A&A, 
425, 577 

Stauber, P., J0rgensen, J. K., van Dishoeck, E. F., Doty, S. D., & Benz, A. O. 2006, A&A, 
453, 555 

Sturm, B., et al. 2010, A&A, 518, L129 

Szymczak, M., Pillai, T., & Menten, K. M. 2005, A&A, 434, 613 

Tafalla, M., Mardones, D., Myers, P. C, CaseUi, P., Bachiller, R., & Benson, P. J. 1998, 
ApJ, 504, 900 

Tafalla, M., Santiago-Garcia, J., Myers, P. C, CaseUi, P., Walmsley, C. M., & Crapsi, A. 
2006, A&A, 455, 577 

Takahashi, T., Silk, J., & HoUenbach, D. J. 1983, ApJ, 275, 145 

Tamura, M., Gatley, I., Waller, W., & Werner, M. W. 1991, ApJ, 374, L25 

Thi, W. F., et al. 2001, ApJ, 561, 1074 

Tielens, A. G. G. M., & Hagen, W. 1982, A&A, 114, 245 



Valiron, P., Wernli, M., Faure, A., Wiesenfeld, L., Rist, C, Kedzuch, S., & Noga, J. 2008, 
J. Chem. Phys., 129 

van der Tak, F. F. S., Black, J. H., Schoier, F. L., Jansen, D. J., & van Dishoeck, E. F. 2007, 
A&A, 468, 627 

van der Tak, F. F. S., van Dishoeck, E. F., & Caselli, P. 2000a, A&A, 361, 327 

van der Tak, F. F. S., van Dishoeck, E. F., Evans, II, N. J., Bakker, E. J., & Blake, G. A. 
1999, ApJ, 522, 991 

van der Tak, F. F. S., van Dishoeck, E. F., Evans, II, N. J., & Blake, G. A. 2000b, ApJ, 537, 
283 

van der Tak, F. F. S., Walmsley, C. M., Herpin, F., & CeccareUi, C. 2006, A&A, 447, 1011 
van der Tak, F. F. S., et al. 2010, A&A, 518, L107 
van der Werf, P. P., et al. 2010, A&A, 518, L42 
van Dishoeck, E. F. 2004, ARA&A, 42, 119 

van Dishoeck, E. F., Blake, G. A., Jansen, D. J., & Groesbeck, T. D. 1995, ApJ, 447, 760 
van Dishoeck, E. F., & Helmich, F. P. 1996, A&A, 315, L177 

van Kempen, T. A., Doty, S. D., van Dishoeck, E. F., Hogerheijde, M. R., & J0rgensen, 
J. K. 2008, A&A, 487, 975 

van Kempen, T. A., van Dishoeck, E. F., Hogerheijde, M. R., & Giisten, R. 2009a, A&A, 
508, 259 

van Kempen, T. A., van Dishoeck, E. F., Salter, D. M., Hogerheijde, M. R., J0rgensen, J. K., 
& Boogert, A. C. A. 2009b, A&A, 498, 167 

van Kempen, T. A., Wilner, D., & Gurwell, M. 2009c, ApJ, 706, L22 

van Kempen, T. A., et al. 2009d, A&A, 501, 633 

— . 2009e, A&A, 507, 1425 

— . 2010a, A&A, 518, L128 

— . 2010b, A&A, 518, L121 



-70- 

van Zadelhoff, G., et al. 2002, A&A, 395, 373 

Vejby-Christensen, L., Andersen, L. H., Heber, O., Kella, D., Pedersen, H. B., Schmidt, 
H. T., & Zajfman, D. 1997, ApJ, 483, 531 

Velusamy, T., Langer, W. D., & Marsh, K. A. 2007, ApJ, 668, L159 

Visser, R., , & et al. 2011, A&A in prep. 

Visser, R., van Dishoeck, E. F., Doty, S. D., & DuUemond, C. R 2009, A&A, 495, 881 

Viti, S., Roueff, E., Hartquist, T. W., Pineau des Forets, G., & WiUiams, D. A. 2001, A&A, 
370, 557 

Viti, S., & Wihiams, D. A. 1999, MNRAS, 305, 755 

Walmsley, M., & van der Tak, F. 2005, in ESA Special Pubhcation, Vol. 577, ESA Special 
Publication, ed. A. Wilson, 55-60 

Wampfier, S. F., et al. 2010, A&A, 521, L36 

Watson, D. M., et al. 2007, Nature, 448, 1026 

Westley, M. S., Baragiola, R. A., Johnson, R. E., & Baratta, G. A. 1995, Nature, 373, 405 
Whittet, D. C. B., ed. 2003, Dust in the galactic environment 

Whittet, D. C. B., Bode, M. F., Longmorc, A. J., Adamson, A. J., McFadzean, A. D., 
Aitken, D. K., & Roche, P. F. 1988, MNRAS, 233, 321 

Wiesenfeld, L., & Faure, A. 2010, Phys. Rev. A, 82 

Wilking, B. A., Greene, T. P., Lada, C. J., Meyer, M. R., & Young, E. T. 1992, ApJ, 397, 
520 

Wilson, B. A., Dame, T. M., Masheder, M. R. W., & Thaddeus, P. 2005, A&A, 430, 523 
Wood, D. O. S., & Churchwell, E. 1989, ApJ, 340, 265 

Woodall, J., Agiindez, M., Markwick-Kemper, A. J., & Millar, T. J. 2007, A&A, 466, 1197 

Wyrowski, F., Mcntcn, K. M., Glisten, R., & Belloche, A. 2010a, A&A, 518, A26 

Wyrowski, F., Menten, K. M., Schilke, P., Thorwirth, S., Giisten, R., & Bergman, P. 2006, 
A&A, 454, L91 



-71 - 

Wyrowski, F., et al. 2010b, A&A, 521, L34 
Yang, C.-H., et al. 2010, J. Chem. Phys., 11, 11 
Yildiz, U. A., et al. 2010, A&A, 521, L40 



This preprint was prepared with the A AS macros v5.2. 



